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2.2 Effect of Reduced Direct Solar Illumination 



R=l-D(l-Spa)/A (4] 

that is independent of both location in the scene and the spectral properties of the surface. R is of course dependent 
on the atmospheric properties, particularly the aerosol amount and type, as well as the solar and viewing geometries. 

A typical calculation of R using MODTRAN appears in Figure 2. Since the Equation (4) result is close to 
the simpler expression 1-D/A that neglects Spa, the replacement of pe with the scene-average pa should be a 
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the minimum distance or angle, or equivalently the maximum angle cosine. 
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Figure 4. Road detection results. Left-hand images are the angle cosine outputs; right-hand images are the results 
after thresholding. 

The classification results are shown in Figures 6a-d. Enlargements of Figures 5 and 6a-d that focus o 
the main cluster of buildings appear in Figures 7 and 8a-d. The color scheme in the classification maps is green = 
grass, ochre = dirt, pink = roof, dark blue = pond, black = panel, gray-blue = road/driveway. 



Figure 6a-d. ARM site classification maps. 
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Results from the Vicarious Calibration of ADEOS/AVNIR and the Visible and 



confirm the validity of the difference, vicarious calibration based on the radiance-based method was 
conducted using Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data. Also, an effect due to 
the spectral refi*active index of aerosol on the estimation of the ToA radiance was confirmed, together 
with smear-like image defect in the AVNIR imagery data. 
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other error sources, in particular, bidirectional reflectance distribution function (BRDF) effect, 
smear-like image defect observed in the ADEOS/AVNIR is described. 



A Brief Description of Previous Vicarious Calibration Results Based on the Reflectance-Based 
Method 
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3. Vicarious Calibration Using the Radiance-Based Method 
3.1 AVIRIS Data Used 


Figure 1 shows the AVIRIS-derived AVNIR and the real AVNIR data. The instantaneous field of 
view (IFOV) and spectral response are adjusted to the AVNIR. A smear-like image defect is seen in all the 
bands of the real AVNIR images. 
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accurately estimate molecule optical depth (Rayleigh Scattering), optical depth of ozone, and water vapor 
compared to that for aerosol optical depth (Mie Scattering). One of the largest error sources is estimation 
of the refractive index of aerosol, which depends on aerosol type and size distribution of aerosol particles. 







Figure 2. Spectral Characteristics of the Real Part of the Refractive Index. 
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(AVIRIS altitude) (bottom part of the table). 



3.4 Comparison of ToA Radiance 
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Figure 5. Transparency of the atmosphere from the 20 km altitude (AVIRIS altitude) to the top of the atmosphere. 
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coincidence so that it may concluded that AVIRIS provides a qualified radiance. Moreover, smear-like 
image defect observed in the real AVNIR imagery data is significant. 
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SPECTRAL UNMIXING OF VEGETATION, SOIL AND DRY CARBON IN ARID REGIONS: 
COMPARING MULTISPECTRAL AND HYPERSPECTRAL OBSERVATIONS 
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where p(X)e is the reflectance of each land-cover endmember (e) at wavelength X. The sub-pixel cover fraction (Ce) 
of each land-cover endmember may be PV, NPV, bare soil or other constituents. Solving for the sub-pixel cover 
fractions (C^) therefore requires that the observations (in this case, reflectance or p(X)pi„ei) contain enough 
information to solve a set of linear equations, each of the form of equation ( 1 ) but at a different wavelength (A,). 
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area, depending on the spatial and temporal variability of vegetation, soils, rocks and other features. Equally 
important, vegetation spectral properties collected at the leaf or branch level, or very close to the top of a plant 
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derive uncertainty estimates of the sub-pixel cover fraction values (Figure 2). The model, AutoMCU, is based on a 
code first developed only for the shortwave-IR spectral region (AutoSWlR; Asner and Lobell 2000), but it is more 
general in that any combination of optical wavelengths can be used in the unmixing process. AutoMCU uses three 
spectral endmember “bundles”, derived from field measurements, to decompose each image pixel using equation 
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15-30m and -20m, respectively. All of the spectral wavelength sampling and signature permutations compared in 
this study are presented in 
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presence-absence as well as in our visual surveys of the area (Figure 1). NPV was extremely high at the grassland 
site (mean = 36.7%, s.d. = 13.3%), with values ranging from 10-70% along the 300m transect. These estimates 
concur with recent studies conducted in the same or similar areas (Schlesinger et al. 1996, Asner et al. 2000). 

Field Spectroscopy 

Field spectral data collections yielded , 


o 

d 






If! 


- o 

t« • 
rri 


ed 


dauepayan pan 



^ s ^ s ^ s g 

o o o d o d d 

aouapaya^ 



S| 

•<? 


=i| 


§.§1 
1^1 
J s o ^ 

J &> 'g)-C .S 
12 

I "S) g 

1 S 

g 


o 

ill 

6 g 

(D 

i’f 


XJ 

< 'B 


aauo)aaya}{ 


aou«)oaya;i 


g d 

a p 

2 d> 

S B 

lx 

.a|| 

(X -S a. 




_ X3 

o> B 
’B IX 

Sp 


t/5 

§"1 

CO 



d 

'S Sx 

P (30 

On ^ 
w On 

(N 2 

>> 
•a 

" T3 

§ § 

^ .b 

4-4 ^ 

6 O & 
S' W o5 § 

S)-e 


CA 

cO 

O 

.1 

TD 

(U 

.3 

1 

O 

a> 

d> 

6 


B .2 

^ ■■ 


o> 
<N 

S'? 
§ g 

2 PQ 

I" 

.s 


1 & 
s 

e-’S g.> 

I a 

^ o-s a 
I . 

o 
o 

CA 


'2|'^ 


(U 

T3 




(D 


(L> 


g 

& 


O 

CA 

a 


5^ 

^ ^ ^ 
2 CA a 
H 2 .2 

^ 3 bo 
• a <l> 

(/3 S3 »-f 

cA X4 rs 

^ O C/3 
2 <D (D 

I 

|:sS 


o 

a 2 ^ 

s 


a 

(/) 

CT 
»n 
(D <u 

5^0 £ 
2 ^ 


l.l 

8 fe 


CA 

1 

O 

o 


I 

g 

§ 

K 

CJ 

i 

I 

i 


s 

2 Ia 

^ d 

;g ^ 

u ^ 

p 

• _r 
CO a 

VO o 

§,1 
E c 
.3 

JJ 

3 


2 


I 

I 


1 


> 

I 

s 

% 

0 

1 
'I 

o 

S 

O 

(D 

Oh 


4-4 

<U ® , — . 

<D 4h 4) T3 
w fl d) 

<D B C 
-8 'S g . 

T3 p4 

0) o> ^ ^ 
CJ ^ Q 

1 I O xo 
tfci -8 ^ 

03 O 

2 ^ 

^ lo S tx 

^ S &'3T 

8 -2 -S ,S 

g o 2 

o § 

cS.§ fe| 
ai3 I § 

‘S 8 I &> 

S "2 

« ^ 8 


C/3 

2 

c 

(D 

rS 

, .2 
</5 d) 

on “ 
2^ 

5 § 


d> 

XJ 

H 


d) 


. 

3 S' c3 'g 
‘gx w)S J 

|sl-| 

^ u “ 

p ^ ^ ^ 
2 ^ 
00 X CJ 
a ^ o <D 

(50 g ^ ^ 

2 Y o 

^ w bX)-“ 

-*-j c/a d) 4) 

bbS p p 

l» 2 f 

a « § a 

^73 ai -2 
^ >x§ g> 

® I dS 

e8 3 B 

fsl p M « 
►Y r;3 d) <D 

3 

1-a S I 

s S a § 

:i o5 YJ o 


VO 

CN| 







0 1 o.taU i n in t. . . 1 

0.6 1.0 1.6 2.0 2.6 0.6 1.0 1.6 2.0 2.5 0.5 1.0 1.5 2.0 2.6 

Wavelength (ftm} Wavelength (/xm) Wavelength 0*m) 

Figure 7 . AVIRIS data of shrubland (upper panels) and grassland (lower panels) study areas convolved to satellite sensor 
channels: (A) Landat TM, (B) MODIS, (C) ASTER, (D) Landsat TM, (E) MODIS, (F) ASTER 
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work presented by Asner and Lobell (2000), and indicate that the tied-SWIR2 spectra are a means for estimating the 
three dominant surface constituents in two dramatically different arid ecosystems. 



appropriate for linear spectral mixture analysis. 
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7. Conclusions 
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1.3 [ im ) to decompose image pixels into PV, NPV and bare soil covers in an arid shrubland and grassland ecosystem 
(Table 2). Shortwave-IR2 (SWIR2) measurements in the 2.0-2.3 jim range showed the distinctive differences in PV, 
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APPLICATION OF RADIATIVE TRANSFER THEORY TO 
ATMOSPHERIC CORRECTION OF AVIRIS Data 
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3. Bio-Optical Model Development 

To develop analytical algorithms for these case 2 waters it is necessary to develop an appropriate bio- 
optical model linking the water constituents to the inherent optical properties (lOP) and to link the lOP to the 
subsurface irradiance reflectance. In this study, we used the linear backscattering albedo model (Gordon et. al, 
1975), which is being considered the most appropriate model for case 2 waters (Dekker et. al. 1994). 














































o 


The absorption and scattering coefficients of the main water quality parameters collected for all stations 
were averaged to produce a set of representative values (Figure 2) and were used as input for calculating R(O-) for 
different dates. No significant changes in concentrations of water quality parameters were observed in the field data 
collected in 1999 and 2000. The shapes of the absorption curves reflect the composition and concentration of 
different optical water quality parameters (i.e., phytoplankton). 



Using equation 1, the modeled R(0-) were generated to represent the spectra for the range of chlorophyll-a 
concentrations calculated in the laboratory for the individual station sampled on May 15 ’00. Figure 3 depicts the 
modeled R(0-) versus measured R(0-) for station 2 (Atlantic Highlands) sampled simultaneously with the AVIRIS 
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RCGb INDEX: A TOOL FOR MAPPING THE DEGREE OF WEATHERING 
IN THE TROPICAL SOILS IN BRAZIL 



The hyperspectral sensor now in operation is AVIRIS (Airborne Visible/InfraRed Imaging Spectrometer), 
property of the Jet Propulsion Laboratory (JPL/NASA). It is an airborne system composed of four spectrometers that 
give images continually scan the reflected optical spectrum from 0.4 to 2.5 |xm, in 224 bands, with widths of 10 nm 
approximately. 



three bands ratio (Gao et al., 1993). The method uses the spectra of pure minerals of the JPL spectral library (Grove 
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spectral math, using 50% of e ach mineral for the formation of a spectrum, as it can be s een in Figme 3. 


Figure 3. Spectra Kaolinite Well Ordered PS-1 A (Kaolinite), Gibbsite Synthetic OH-3A (Gibbsite) (Grove et al, 
1992) and the integration of 50% of each one to generate the spectrum Kaol+Gibb. 



ro 




several wavelengths, adapted to AVIRIS bands. 



Where bl, b2, b3, b4, b5 and b6 correspond respectively to the reflectance values in the wavelengths of 
2,127 }xm (band 186), of 2,176 pm (band 191), of 2,196 pm (band 193), of 2,226 pm (band 196), of 2,266 pm (band 
200) and of 2,286 pm (band 202) (in the case AVIRIS data) or 2,335 pm (equivalent the band 207) (if laboratory 
spectra are used). The use of three band ratio procedure, besides the best Y scale adjustment, minimizes the noise 
found in the AVIRIS spectra. 


Figure 4 presents the two integrated spectra, with the continuum removed, showing the six bands used in 
the RCGb spectral index. 
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Of the 56 spectra, we removed four that presented opaque minerals. The influence of those minerals in the 
determination of the RCGb spectral index was discussed by Baptista et al. (2000a), using samples of the municipal 




district of Niquelandia, Goias. The 52 diffuse reflectance spectra obtained in laboratory were resampled to the 



mgn correlation coemcient as it can be seen in Figure 5. 

Ikaol allows only estimating the kaolinite contents in relation to the fraction clay in a satisfactory way. The 
inclusion and the analysis of another parameters can drive a more precise estimate method for the kaolinite content 
of the latosols, by means of r eflectance spectra (Madeira Netto, 1993). 
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This work showed that the RCGb index could minimize the time expended in the analysis and in the 
recognition process and mapping soils. To perform those two study areas in the traditional way it would be 
necessary to adopt a sampling matrix with regular spacing covering the whole area, which would increase the c 
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In the intimate mixtures series among the goethite and the pimelite accomplished in a laboratory 
a 10% intervals abundance was also used (Figure 4). The SCM in the pimelite absence also 
presents high values demonstrating a high correlation among the spectra. The coefficients have 
linear behavior (Figure 5, 6, 7, and 8). 



4. Abundance of Garnerite Analysis when Mixed with Asbolan 

The gameritic facie of the mine has asbolan veins. In spite of the presence of the asbolan being 
restricted to a specific area, the use of that mixture is interesting to demonstrate another type of 
spectral mixture behavior. 
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Wavelength (pm) Wavelength (pm) 

Figure 1-Simulation of the macroscopic mixtures between the pimelite and the goethite with 
abundance interval of 10%: a) without normalization (removal of the continuous) and b) with 
normalization. 
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Figure 4-Simulation of intimate mixture between the pimelite and the goethite with interval 
abundance among the curves of 10%; a) reflectance (offset for clarity) and b) continuum rem 
reflectance. 
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Figure 6-Linear coefficient Figure 7— Angular coefficient behavior for 

behavior for the intimate mixture the intimate mixture between pimelite and 

between pimelite and goethite goethite 
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Figure 10-a) Angular coefficient (B), linear (A) and it scale (E) behavior for the macroscopic 
mixtures simulation between pimelite and asbolan b) SCM among pimelite and the relative 
curves to the macroscopic mixtures between pimelite and goethite 
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Analysis Absorption Band Positioning: A New Method for Hyperspectral 

Image Treatment 
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Figure 1 - Wavelength image of the absorption band. 
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Figure 4 - Spectral curves refer to green vegetation (0,6849 pm) and nonphotosynthetic 
vegetation (0,6754 pm): a) reflectance and b) continuum removed reflectance. 
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Figure 6 -Spectral curve sequence between goethite and hematite (0,46 |o,m and 0,62 




Another interesting application of AABP in the study area is in the characterization of 
the gamerite main spectra representative of the facies with a concentration of nickelliferous 
silicated minerals. The term gamerite has a generic connotation to define a mixture of 
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Figure 8 - AAPB image relative to the Mg-OH band. 
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information can be used to evaluate the degree of correlation among the materials 
analysis. 



Where: 

MSr - Sums of squares mean regression 
MSd - Sum of squares mean deviation 
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These features used the value of degrees of freedom in the denominator equal to one. 

Compared to the manual analysis, it is observed that the estimate of 2.5% was quite 
efficient to identify the areas for both features (Figure 2 and Figure 3). The other 
significance levels present an overestimated area. That characteristic is due to the high 
spectral correlation of the scene. 






Figure 3. Comparison of the areas defined for kaolinite presence, using the SIM method for 

the significance levels 2.5%, 5% and 10%. 
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The method has a similarity parameter that completely eliminates the false positives present 
in the algorithms of SAM, SFF and Tricorder and still maximizes the values of the areas 
where the sought material exists. 
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Figure 5 Scatter plot among SCM x SIM for the features: a) pimelite and b) green 
vegetation. The curves are segmented in agreement with the degrees of significance: 2.5% 
(red), 5% (cyan) and 10%(blue). The areas in yellow have a relative value higher than 10%. 




mployment of the Multiple Endmember Spectral Mixture Analysis 



describes each pixel. A routine was developed in IDL, which accomplishes that 
comparison of the models of mixtures defined for each pixel the most appropriate. This 
method was tested in the Niquelandia (Brazil) mine and showed excellent results being 
in agreement with the field. 

Keywords: Spectral mixture, MESMA, imaging spectroscopy. 
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the use of spectral classifiers, for instance, the Spectral Correlation Mapper (SCM) method 
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a) Pimelite - Saponite b) Goethite 



Figure 2 - The rms image 










Weighting coefficient 
on continuum -removed spectra 

Figure 4 - Relationship of the scale coefficient with the weighting coefficient on continumn- 
removed spectra (normalized absorption band area). 






a) Pimelite - Saponite ^ Goethite 



Figure 6 -Disposition of the models selected by the MESMA considering 

continuum-removed spectra. 
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INVESTIGATION OF IMAGING SPECTROSCOPY FOR DISCRIMINATING URBAN 
LAND COVERS AND SURFACE MATERIALS 
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This preliminary study demonstrated that the imaging spectroscopy approach with low-altitude AVIRIS 
data and field spectra makes it possible to discriminate the urban land covers and surface materials. The low-altitude 
AVIRIS data with 2.9 m spatial resolution successfully differentiates many surface materials in the urban scene. 
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The objective of this work is to discuss the basic problem of NDVI variations considering the atmospheric 
effects and the band positioning and bandwidth influences. For this purpose, the position and width of narrow and 
broad R and NIR bands, including those from some multispectral sensors, were simulated from AVIRIS (Airborne 
Visible/InfraRed Image Spectrometer) data collected in central Brazil, and NDVI values were calculated and 
compared from apparent and surface reflectance images. 



METHODOLOGY 



where is the spectral reflectance, is the spectral radiance, E;^ is the spectral irradiance on the top of the 
atmosphere (complementary information accompanying the radiance-image) and 0s is the solar zenith angle obtained 
to the central point of the scene, considering the date and time of image acquisition. In the present investigation, 0s 
was of 3 1® . 



absorption bands in the blue and red intervals, the reflectance peak in the green, the absorption features due to leaf 
water in 960 and 1200 nm, and the lignin/cellulose features around 2100 and 2300 nm. The S spectrum is from 
Latossolo Vermelho-Escuro (oxisol). Its surface reflectance spectrum displays a broad absorption band at 900 nm. 



due to hematite, and a well-defined narrow feature at 2200 nm due to kaolinite. These GV and S spectral 
characteristics are affected in the visible interval by atmospheric scattering and in the NIR interval by atmospheric 
absorption. The strongest water vapor absorption bands are located at 1400 and 1900 nm, spectral intervals not used 



Figure 1 - AVIRIS-derived surface (continuous line) and apparent reflectance (dashed line) curves of 
green vegetation (riparian forest). The main atmosphere constituents responsible by well-defined 
absorption bands are indicated. The horizontal bars represent the nominal intervals of R and NIR bands 
of the sensors listed on the right side of the illustration. 


Table 1. Spectral nominal intervals of the R and NIR bands of some orbital multispectral 
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the two major intervals of atmospheric absorption. 
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Figure 5 - NDVI difference curves derived from surface (ps) and apparent 
(pa) reflectance spectra of green vegetation (GV) and soil (S), as a function 
of the positioning of the R and MR bands of the sensors listed in Table 1. 
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SPECTRAL CLASSIFICATION OF SIMILAR MATERIALS USING THE TETRACORDER 
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capacity in the watershed therefore requires knowledge of the distribution and abundance of all three minerals. Thus, 
solving the Calcite-Epidote-Chlorite problem became a primary interest for the authors of this study. 
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quantitative unmixing approach based on optical constants (Hapke, 1993; Mustard and Sunshine, 1999); and the 
second is the less elegant empirical approach of comparing spectra of an assortment of fine-grained laboratory 
mixtures. Given the present limitations of modem microprocessors, this study utilizes the latter approach. 
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misidentification becomes much less likely. 
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The next step was application to the AVIRIS Animas Watershed flightlines. These two flight lines were 
acquired on June 18, 1996 at 9:30 AM and cover the region from Durango, Colorado to the headwaters of the 




correlations to the pure endmembers and were therefore not considered in the final analyses. As chlorite and 



muscovite had already been encountered together in a different location (Swayze 1997), a chlorite-muscovite mixture 
was already in the USGS Tetracorder (v3.4a8) reference library (Dalton, 1998). 

Four primarily northwest-southeast trending zones of exposed propylitic alteration are visible in Figure 5. The ridge 
of Hazelton Mountain in the upper right comer is highlighted in a profusion of light and dark blue, yellow-green and 



Figure 4. Georeferenced image of a portion of the Animas River Watershed study area. Tetracorder results 
are shown in color. The town of Silverton, Colorado is near the top center; the Animas River bisects the image 
vertically. This image has been engineering-corrected for roU, pitch, yaw, forward velocity, and smile, as well 
as orthorectified against a digital elevation model to remove topographic distortions. 
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the case of mixed pixels of epidote, calcite and/or chlorite, together with the muscovite signature at 2.2 pm, nothing 
further could be determined about the composition and these are the orange pixels in Figures 4 and 5. 
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Clark, R.N., G.A. Swayze, T.V.V. King, K.E. Live, J.B. Dalton, and R.F. Kokaly, Expert system features 
identification rules for reflectance (and emittance) spectroscopy analysis I: Visible to near-infrared detection of 
minerals, organics, vegetation, water, amorphous and other materials, J. Geophys. Res., 2001. 


:2 

3 ^ 

Oh ^ 

S S 

<D K 


C/5 

<i 

d 

'O 

§ 


o 

(/) 

s . 

^ OO 


C/3 


8 2 


O 

c3 

> 

CO 




fO 

00 

I 

0\ 


u 


3 a 
< 
(D 


t 

Oh 




^ s 'I 
^ ^ ^ 

I .9 ^ 


s 


oi) ,S 


<CJ 

op Co 
.9 ..s* 


c 

2 

> 

> 

H 

PQ 


G 

O 

13 

Q 




‘S s 

K.q 

•a 

2 
(U 
G 
(U 


Q 
S 
C 
K 
W) ^ 


G 

4h 

O 

G 

*G 

G 


<U 

CO 

CO 

< 

d 

§ 

<lT 

o 


J? 

1 

Ci) 

CJ 

I 

1 

u 

w 

o 

U 

o 

(L) 

Gh 

CO 

00 

2 

5 ; 


r< 


> 

C 

>* ^ 
^ Gh 

T3 


>% 

Xi 

T3 

<D 

Xi 


1 : 

u ^ 
Z 2 

2 s 
S 


O 

o 
o 

<N 

g O 
t3 0 


of 

> 

o 

PQ 


X ^ 
Q .G 

bB 

” ‘g 

X! 


> 

> 

H -o 
§ 

PQ M 


7^ Cj 

G Q 

e 
o 


G 

O 

'*-» 

'g 

Q 




2 
<D 
G 
(U 
W) ^ 


W 

2 CO 

>;2 

■ l-H 

> 


G 

2 


X 

Xi 

"O 

(D 


(D 


d> 

Vh O 
(D 
> 


CO 


G 


(L) 


2 -S 


(D 

G 

a> 

Q 

I 

o 

.0 

§ 


(D 

1 

O 

PQ 

U 

Z 

2 

o 


G 


^4H 

O g 

G Cl 
.1 ^ 
c3 ^ 

■S ! 

g Co 
G 

is ^ 

G o 
61^ 
.si 

CO 


1 o 


> 

> 

H 


^ 00 


o 

PQ 


PQ 


'G 

§ 

<D 

O 

, § 

> 

CO 


-St 


G 

o 

4-J 

13 

Q 


"9 

-St 

-SJ 

-G 

CIh 

-G 

Gh 

o 

o 

00 

S 

o 

<L) 

Gh 

CO 

• ^ 
W) 
G 


W) 

G 


I 

(L) 

G 


4h 

O 

(D 

2 

U 


;0 tri 

G 


T3 

§ 


*C 

d) 


G 

S 


§H 

PS 


.3 

o wa 

1 — H 

13 ^ 

.2 G 
bO ^ 

o 'S 

2 S 
2 


.> 

2 

T3 


g J 

u PQ 


St 

§ 

r 

CO 

8 

§> 


S 

^ ■§ 
T3 G 

§ o«a 
M 

^ O 
W) 2 

.pi T3 

c§ a 

K G 
2 T) 
.V G 

JS 

O g 

§ !a 

O G 

® S 

c/a 

Q Gh 


^ O 

o ^ 
3 Ci) 
S Co 

c/a ^ 

T3 


cd 


§ 

S 


CO 


O G ^ 
G PQ 00 

On 


O 

>Ph 


^ X 


On 

On 


G 

O 

c/a 

G 

H 


o 

G 


V, 

O 

o 

(L) 

qS 

Q 


(D 

G 


(D 

o 


(N 

o\ 

S 3 

I <2 

s J 


it 

1 1 

8 -S 

Gh 

00 Co 
T3 

II 

bO S 

•S o 
D :§ 

•-«.* 
o 
s 

St 
St 


§ 

I 

o 

PQ 


T3 


o 0 ^ 


73 

2 

PQ 


X 

(J 

(D 


oa 

O 

S 

S 

CO 


X 

‘C 

ffi 


I 

00 

2oo 

ffl > 2 : 

K/ <C 

r^S 

^ C 

£J 8 

S g 

a 8 


0\ 


00 

VO 


ON 

»o 


fO 

On 

On 


CD 

bD 

73 

•G 

I 

U 


c/a I 


d> 

bO 

73 

•c 

I 

u 

I 

I 

8 

S’ 

0 
St 

§ 

1 

St 

Q 

VJ 

St 

a 

§ 


PQ 

clT 

§■ 

m 



G 

> 

X 

c3 

CD 

n 


K 

.§ 

4 

> bo 


8 ^ 
c< 

G 

O 25 

S 2 


CD 

+-» 

CD 

I 

O 

CD 

Gh 

CO 

bO 

G 

■§> 


73 
CD 
c/a 
G 

s£ 

B 'O 

c^ (D 

^ s 

W rO 


PQ 


> 

CD 

8 

o 

X 


I o 


G rt 

U 4 


s 

t 

. g 
§ ^ 
H? Co 

§ ^ 

CD 73 
X G 

«+-H O 

0 o 

2 c/a 
c« G 
X ■ G 

.9 s 

1 O 

ia S 

o a 
> § 
73 G 

*-~i 

•- S 

e 

^ o 
X tn 

G CD 

-S ^ 

G (D 

PQ 5 

C/3 .G 


2 

2 ^ 
(lT o 2 

^ rvf 

J (D 'V 


73 


X\ 

d 

2 


r. 


CD 

0 2P 


> 


CO 

C/5 

(D G 

C3N 

-a 
m ^ 

< u 

On 

VN 

CO 

CO 

1 

< 

H 

<f ^ 

CON 

0 

2 

g 

CO 

s' 

(lT (P 

^ !g 

Gh 

G 

8 


5:5 

J NO 
73 

§ f* 


s 

€ 2 
CO ^ 

^ -»i.» 

d) Q 

feb 5 

J3 S 

I I' 

•& « I 

hJ pi Co 


<N 

O 



Mast, M.A., P.L. Verplanck, D.B. Yager, W.G. Wright, and DJ. Bove, Natural sources of metals to surface waters in 
the Upper Animas River Watershed, Colorado, Proceedings of the 5th International Conference on Acid Rock 



m 

o 




ASSESSING THE SPATIAL DISTRIBUTION OF PLANT LITTER 
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The objectives of this research were to (1) determine the spectral reflectance of green plants, plant litters, 
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the soils in the trays were saturated with water, allowed to drain freely for 2 h, and weighed before the spectral 



























reflectance was measured again. Relative water content (RWC) was calculated as the water content divided by the 
maximum water content of each samples (i.e., RWC = 0.0 is oven-dry and RWC = 1.0 is saturated). 
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End-members for linear spectral unmixing were selected from the corrected AVIRIS scene. The regions of 
interest (ROI) included: tarmac for asphalt, recently plowed fields with no residue for soil, several ponds for water, 
several roof tops, and fields of winter wheat and hairy vetch for bright green vegetation. Furthermore, we added the 
reflectance spectrum of a shaded Spectralon panel for the shade end-member, and the laboratory reflectance of com 
residue for the residue end-member. The linear spectral unmixing routine in ENVI 3.2 was used with the seven end- 
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Many chemical bonds in plants have vibrational absoq^tions in the 1500 nm to 2500 nm waveleneth region 
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Table 3. Mean NDVI, CAI, and end-member relative weights from the linear spectral unmixing analysis for regions 
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Figure 10. Relative abundance of soil (upper left), residue (upper right), and green vegetation (lower left) end- 
members. NDVI image (lower right). In each image, white represents high values and black represents low values. 
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CALIBRATION AND VEGETATION FIELD SPECTRA COLLECTION 
FOR THE 2000 AVIRIS HAWAII DEPLOYMENT 


o 


a>. 


s g ® 
■£ 

c« is ^ 

“ -2 2 
> S C3 
O S 

' o w 2 

1 4 h « 

® § 3 

w ^ o 
g “ •§ 

•S >n P 


^ c/a XJ 

« 2 > 

•s’^1 

O ^ 

; ^ si 

\ ^ S 


^ (D Q. 
1>3 V-. ^ 

^ ^ c 

IS ^ 

Pu . 
P (/) 

D <D 

(D +J ‘5 

tjD O (D 


B 13 8 ) 

O O (D 

a w) > 


<D o « 
^ <N «i 


Vh O ^ 
^ ^ 
t> 2 ^ 

g §1 

o> o 

Q, *0 o 
a> fO 'S 
f=5 U 

w - 

tLi ^ jrj 
w- > (D 

(D 

2 I PQ 

^38 


8 § 

CL ^ ^ 

^ 'tl* ^ 

, eg c:> 

cJ w 

o <u ^ 

li^. 


CL 2 ’ 

^ S :5 * 

y O c« 

g 0> o 

^ Ph CL 

w) ^ rt7 
p ^ 8 
*53 ® P 

JD p eg 
"P '§ 

.g § 

8 S * 


2^ S 

S CL T 3 

I O ^ 

jz; (u o) , 
Q w P 

w >,w 
t /5 e« »irt 

^ c ^ 
-9 -P P 


|aa 

■•3'Mi 

eg bo c /3 


0 S 13 S?n 
*S > £? 

<D - g c3 

3 I febe 
^0^0 

1 O 1 

eg 

c« 00 ^ 

S 2^ 3 9 
*5 

1 rfc'S 
®'2|'2 

t- Oi (U 

52 w ^ ^ 

) u ^ <u P 

§ <2 ^ i 

c ''^ ^ eg 

0 (D 

't g go a 

1 -I s 2 

ia |-S ^ 

'H ;i3 

p <1^ 9 

2 G *5 

« ^ I 2 

B “ 3 § 

*- iq •§ T 3 

I 2^ ^ 

tJ 22 u ’C 

| 3 Sl 

i>^ ; 

li^.s 

^ rP i 'P 

Ib |s 

^1 = I 

9 3 ‘9 o 
1 1 2 ^ 
s S ^ 2 

^ o o ^ 
^ ^ 3 ^ 

43 <D S i 2 
^ p CP P 

hil <u 8 EP 
ffi .g g I 

a “ -q 

O rL rti *G 


> o C 

3 ® • 
•S ^ c 

o^-d, 

§ II 


fe P Ip 
■•2 3 2 


I' l« 

eg e« o 

f S 5 

bO P 
»p 'T3 ^ 

wb c 


ftH -2 
^ . p 

^ d rS 
X S c« 
di — 

«a CJ _c3 


3 -P H 


+i > <D 

s a S3 

CL CL,P 


eg 53 ^ 
(D ^ *5 


w «N 

’Gj bo ^ ^ 

^ b o c3 
i 2 S S 5 


§ 3 -2 ** 


S ® ^ ffl 

2|«i 

o> 8 S 

h Ph p *p 

g I li 


p I ^ < 

sl |1 

2 2 a ! 

u § c« ! 
H S M ^ 

,-§ 1 : 
a § ^ 

j2 2^ ! 

® i| ! 

1 ^ h < 

P pS O 4 

1b 

I “‘ 5 -^ 

g-® i j 

s ^ 

i^-i ! 

g,r- eg ( 

' " 3 'd § ! 
s « ^ 

I § I i 

G-S g P 
ca *n S) , 
*-• rq P •! 

<2 S S ' 
P 3* p p 
^ ^ ' 
x> P 

IS 

•q -5 ^ 1 

3 > 2 1 

■g -q ' 

I «■ 2 -1 
iJ i 1 ) 1 

C 4 ^ 

^ gi a . 

a -c « 1 
S 5 .ti 

2 ' *53 

iS S ^ ^ 
^ -S (D 
li II e^ 

^ ^ i 

£ 2 ' 
bO-h 'p 

. -c -S ^ . 

I 4:3 eg PC 


2 

eg 

:a 

T3 

g 

u 

(D 

>y 

s 


eg 

’t-i* 

13 

U 

*S 

bO 

(D 

0 

T3 

> 

bB 

g 

2' 

’o> 



2 


>-H 

eg 


bj> 

CL 

«J 



ill 

^ ffi oi 
^ 'S :S 

si I 

.s ^ bS 

jG bfl W 

^ VP L-i 

•p PQ o 



was taken. Papohaku Beach was the widest beach encountered in the islands, and is considered the better of the two 
calibration targets for Molokai. Comparison of spectra from Kalaupapa and Papohaku Beach reveal differences in 
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vegetation mapping. 



Kalaupapa 



Figure 2. Matched filter results for three species found on the Kalaupapa Peninsula and 

topside Molokai. 




Table 1. Calibration Targets 







Typically, the VRAT is planted in a ridge till, monoculture com and is watered by a central pivot irrigation 
system on a three-day period. Since water stress can increase spectral reflectance from com leaves (Wooley, 1971), 
it is important that the N-application plots be adequately watered so that only nutrient-related stress will 
predominate. Figure 1 shows imagery of the USD A VRAT site with the fertilizer amounts for each block shown. 



work was supported by the NASA Earth Science Applications Directorate under contract NAS 13-650 at the John C. 
Stennis Space Center, Mississippi. 



Low-altitude AVIRIS hyperspectral imagery was acquired over the Shelton, Nebraska, VRAT site on 
July 22, 1999. The overflight produced 3-meter pixels with 224 spectral bands. Ground personnel supported the 



propagation, neural network for training. 



100-kg/ha block on the eastern margin of the field visually shows stress. This stress, however, can be somewhat 
ascribed to a misapplication of the fertilizer - but not all. 



The AVIRIS NDVI image also picks up a stressed area of vegetation, seen as a darker-toned band on the 
north-central periphery of the field. The irrigator is turned off over this area to create water stress effects in the con 
This deliberately stressed region is used in this study to assess whether a given algorithm, or neural network result, 
can distinguish nutrient stress from water stress effects in com. 
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Figure 3. NDVI Algorithm and ANN Data Plot 


The graph on the left in Figure 3 shows the mean NDVI values plotted against the N-treatment blocks from 
which the NDVI values were extracted. Error bars are also plotted and represent an excursion of one standard 
deviation above and below the mean value. In image terms, a given mean NDVI value is really a mean brightness 
associated with an N-treatment plot. The error bars, then, express the dominant range of brightness values that is 
associated with an N-treatment plot. 
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4. SUMMARY 
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Farr View Consulting/Space Science Institute, Boulder, Colorado (farrand@colorado.edu) 
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being the QR factorization rather than SVD. 
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where: R’ = the inverse of the sample correlation matrix of the hyperspectral scene. Since R is ill-conditioned 

for hyperspectral images, it is difficult to calculate R ^ To circumvent this difficulty, the sample correlation matrix 
can be approximated by: 
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►tained from the Clark et al. (1990) band mapping algorithm. In using either approach, the analyst must exercise 


care to use only those channels covering the absorption band of interest. In an early trial, the author used most of 
the AVIRIS spectrometer D channels with the ENVI “spectral feature fitting” algorithm and the results were 




le, G. and A.F.H. Goetz, 1993, “Terrestrial imaging spectrometry: Current status, future trends”. Remote Sens. 
vol. 44, pp. 117-126. 




MEASUREMENT OF COLOR IN HYPERSPECTRAL IMAGES (AVIRIS) 
USING THE CIE (COMMISSION INTERNATIONALE D’ECLAIRAGE) 







physiognomy of the vegetation. However, fires alone do not explain the abrupt limit among the 
gabroic rocks (mafic) and the ultramafic substrate. 
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And the chromaticity coordinate x, y e z is expressed by: 



method proposed by Green (1990) was specifically developed for AVIRIS hyperspectral images. 




The Green method provides an estimate of the atmospheric parameters and a calculation of the 
apparent reflectance of the surface using the code of radioactive transfer together with a model of 
non-linear adjustment for the square mean (Green, 1991, Green et al, 1993a,b). 


« o o 

« « ^ g 


c/) 


4h 

O 

<D 

a 


cd 


S 8 ^ ^ 

&, 'o o 

Sj ^ S 

S ^ 

O t/j 
> (D ^ 

c 3 O O 

§ O 

•Sf> g a « 

^ (/3 

.a M 2 ^ 

^ ^ § I 

(1-. 2 o 

^ ^ o 


M 


§ 


!l> 


.9 ^ 

" •»-( 


{=! 

O 

• 

o 

o 

o 

*a 

CT) 

O 


;h 

(D 

cK 

c 


N 


s S ^ 

cd ^ (U « 
(D S3 cd 

n:^ ^ B 

S .3 S 

^ O o ^ 

>> Cl( (D rH 

^ S ^ ‘42 
C/D cd > cd 

I i a & 
I'i 

^ (D S 

^ ^ 

.O ^+H 

« 2 ° 
T3 (D O g 

S 2 « 8 

cd cd li <D 
N > ^ §3 


O 

a, 

(/) 

. 

a ^ 
o 

a ^ 

o ■S 

>H O 
4 ) 

II 

(VI 

H 8^ 


^ a 

li^l 

§ Oh « OJ 
M (-H (/) ^ 

m § ^ i 

W • rj O P 

n td P 
^ ts ^ P 

o bX)^ c« 
^ U C3 

g ^ hX 


o ^ 

■« a 
M ^ 


■c 

cd 
"d 

(/) *S 
P 

Bn'S 

a ^ 

Id 


O 

r-t 

5 t/) 

O (D 

'42 a 

cd ' 
o ^ 

^ *a 

I «’ 

rs 


C/D 

<D 

C/D 

P 

td 


<D 


Cd 


C/D 

O >-i ‘i 
O cj i 

ll 1 

o ' 

O ^ 


d S 
E ^ 
§ 


a. 

o 

0 

1 

I, 


p 


0 

d <+H 
cd ^ 

S ^ 

1 

»—H -*-» 

^ bo k 
« o' 
<u > 

2 ^ 


*0 
o 

PQ 

(D 

- t3 

S 3 (D 

1“ & 
2 o 

MH c/i 










Kfl I 

§3 s 

M CS .2 

'C « 

S .2 § 


For the Niquelandia area it was possible to separate for the parameters of CIE vegetation 
patterns in soils originated from ultramafic or mafic rocks. 
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Rejection Basin Detection Using the Integral Spectral Analysis (ISA) Method, 

Niquelandia, Brazil 



In 1973 the Niquel Tocantins Company built a metallurgic complex to produce 
nickel carbonate from lateritic ore. Nickel carbonate is transported to Sao Paulo, where 
metallic nickel is produced. The wastes are stockpiled in sedimentation basins around the 
company. Mainly opaque minerals (maghemite and magnetite) constitute the wastes. 

The AVERIS image shows the mill, rejection basins and the dam (Figure 1). The 





terms of image quality (Green et al, 1988.). 





Figure 3 - MNF image to the study area. 
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Our research to date has focused on two primary objectives: development of a regionally specific spectral library for 
urban areas, and using advanced techniques to map urban materials from hyperspectral AVIRIS da ta 
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range of pixels depending on the size of the target. Once extracted from the image each spectrum was given a unique 
identifier in the spectral library. 
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MESMA will also be the technique used in road condition assessment, although our initial emphasis in this project is 
on road extraction. Once the road extraction technique is completed, a spectral library including roads of various 
conditions will be utilized within MESMA to map road quality. 





library is limited in that endmembers were image derived-homogeneous pixels within the image were averaged to 
build the endmember library. Due to the 3.9-meter spatial resolution of the imagery, these endmembers inherently 



contain mixtures of different surface materials. Purer endmembers may be obtained by using a field spectrometer to 
collect known surface types. Therefore, the next major step in this project will be to develop a reference spectral 
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Spectral Library Development 



Figure 1. Methodology for creating the urban spectral library. Homogeneous areas were first identified in the 
field and then identified in the AVIRIS image to be included in the library. 






Pan-chromatic image of urban Q-Tree filtered white lines, 

intersection. 

Figure 2. Q-tree filter concept. This approach incorporates spatial information into the procedure to alleviate 
spectral confusion, particularly between composite shingles and roads. 
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OCEANOGRAPHIC AND ATMOSPHERIC RETRIEVALS FROM AVIRIS 
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. Modulation Transfer Function 





Incosd 


Eq solar irradiance at top of the atmosphere. 

A electromagnetic wavelength. 

V C siirface slope required for specular reflection. 
T two-way atmospheric transmittance. 

psp Fresnel reflectivity of water surface. 



the surface reflectivity is that corresponding to the air-seawater interface. When environmental 





Figure 1. Raw data image of band 54 excised from scene 4 of AVIRIS flight f970410t01p02r02. 
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We fit the total-slope variance, transmittance, and path radiance in the above equation to 
the ensemble. The data and the goodness of the fit for the April 10“’ case at several wavelengths 
are shown in Figure 2. The fitted values for the transmittance and wave slope variance are shown 
as a fimction of wavelength in Figure 3. In this case the total wave slope variance is 0.06. Using 




buoy data to estimate that portion of the total wave-slope variance produced by the large waves, 
we can estimate the variance of the small-scale wave slopes. 
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Figure 5. Buoy and optically-derived wave spectra from AVIRIS data acquired on April lO*. 






Surface reflectance of direct solar radiation also has been computed using the infrared 
bands. The physical parameters such as wave slope computed for these data are, of course, 
electromagnetic wavelength independent. They can be used to remove the effects of glint in the 
visible regime with the Fresnel reflectivity appropriate for the wavelength. 




We have removed atmospheric and solar glinting effects to estimate the upwelling water- 
column radiance or the remote sensing reflectance. This reflectance can be used to to calculate 



The environmental parameters were obtained using the large field of view of AVIRIS for 
homogenous atmospheric and oceanic conditions. If the surface-slope parameters are obtained 
from other sources (e.g., buoy measurements or models), then the atmospheric parameters can be 
retrieved from Equations 1 and 2. This technique enables a local measurement of atmospheric 
parameters that may be generally more applicable, especially for spacebome observations. 


This work was funded by the Internal Research and Development Program of XonTech’s 
Special Studies Division. Buoy data were obtained from the National Data Buoy Center’s web 
site and the AVIRIS data were provided by JPL. 
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HYPERSPECTRAL IMAGE ANALYSIS OF CORAL REEFS IN THE HAWAIIAN ISLANDS 
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Coral bleaching is perhaps the most obvious visual i; 


5 (U 


^ > Cj ^ 

> (U 


Q-» E iil fi Cd 

-2 aj C .2 a e fe 

B ^P!o 

? O O cd ° 


^o=el3 2 


^ ■*“' (/3 

•r c« 2? 
tt a -2 

« g jj 
G O 


rnr;rir'{/aOQ.7^>-i 

|Sh5s«s|8 
^ M 2-^2 ^ 


’S^2g*;;-cdg 

«i§||:ll|i 

ox f=i oS’c« 5-^ 

il.s r 


•o§<Ss-22S’H 

W)o cd 

O .. p '>5 q ^ ^ 


(D C 'lU ^ 

o c2 8 

c+H O G 

0 O ^ 

a o cd 

■ct'^o 

1 I •« “ 

.2 s W 

/Cd Kn rn 


.22 -S 

? u -R "2 


O T3 - 

U ^ g 

s ^ 

S 2 /p 


o ^ 2 

_ D (/) <D 


^c2 g?sg'^^ 

o 2 p S xJ 

o 5 a S “ “ 
o P o > ,S .S kR R 


5§ S " 

“ is 

S S .S S 
S i « P 

^ w D s 
G - S 
- bfi O P 
tii G C .2 


- &‘og a c.g 

-CmOuuS'C(u> 

^C.sr.r-HS^>GD 

oS'tis'S-goS’g 
•l^l'g-^ SS 2 S# 

w.Sc>(D'<^5^2S 

^^3<ug^gjgx3 

Sg’S.sgo^irs 

'P'S® 2-5^2^ ^ 

g s B S ® 1-2 


G 3 JD fc 

^ . o • 

s ^ ^ cB 

Zh .3 (U 

2 I S 

c3 S 3 o § 

U o S g O 

. O (D ^ M 


fc 3 eg I § 

• G G c/5 

t— ( p rj P c/5 

B 3 § J3 2 

^ w) g -g 


sac 

o 23 S 


rj '*“* *^ O ^ 

§ O (/) Cd P 

JZ § g a 

gc/iCd^-- 

O j3 Cm ^ ^ 


g ',£ ‘I 13 3 S ^ S 
3 G S-S^ 

Oc«OL>'POc/5T3(D 


U G 3 

g’l I 

D Gh 

G g 

1 C-2 

g « 33 

«<! s 


health and biodiversity. They also contribute numerous other economic and social benefits, as well as an inherent 
aesthefic appeal. Fortunately, the ecosystems themselves are incredibly resilient and have the ability to recover and 
thrive if adequate management and conservation measures are implemented. However, the complexity of both the 
natural system and the many factors involved does not lend itself to simple solutions, particularly when human 
factors are included in the equation (Hodgson, 1999). 
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that this campaign has made available for coral reef analysis, which includes both significantly human impacted 
ecosystems in the main Hawaiian Islands and the relatively pristine environments of the Northwestern Hawaiian 
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4.2 Kaneohe AVIRIS data 


The flightline used for this project covers the coastal area of the windward shore of Oahu, including the 
entirety of Kaneohe Bay (AVIRIS flightline: fD00412t01 p03_r08). Note that this flightline is predominantly cloud- 
free over the bay (Figure 4). 
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where it can be reflected, absorbed, or transmitted. Reflected energy, however, must pass back through the water 
column, the water-air interface, and the atmosphere before finally reaching the sensor. Ultimately, because the 



shorter wavelengths are strongly affected by Raleigh scattering and wavelengths in the near-infrared and longer are 
significantly absorbed in the water column, the remaining window for analyzing benthic reflectance is limited to 
approximately 400 to 800-nm. 
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Figure 6. Unsupervised classification of Kaneohe Bay AVIRIS data: a) radiance 
image; b) classification output. 





Hochberg, EJ. and M.J. Atkinson, 2000b, “Spectral Reflectance Characteristics of Coral Reef Benthic 
Communities,” Ninth International Coral Reef Symposium, Bali, Indonesia, 23-27 October. 

Hodgson, G 1999, “A Global Assessment of Human Effects on Coral Reefs,” Marine Pollution Bulletin vol 38 
no. 5, pp.:345-355. ’ ‘ ’ 
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by 1-km multispectral satellite 
sensor showing the location of the 
active volcanic vent on the Big 
Island of Hawaii. 






Figure 3. AVIRIS image of the active 
volcanic vent on Hawaii acquired on 
April 26, 2000. 
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Figure 6. AVIRIS image at 2500.9 nm 
showing the emitted radiance of the 
hot lava. 
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Figure 7. AVIRIS spectra spanning the range of major spectral signatures in the April 26, 2000 volcanic 
vent data set. At 2000 nm AVIRIS is saturated by emitted radiance over the hot lava target. 
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Figure 10. Fit of Planck function to AVIRIS spectrum for hot spot A with correction for solar reflected 
radiance. A temperature of 923°K and area fraction of 0.24 was derived. The AVIRIS dynamic range is 
saturated in portions of the spectrum beyond 1500 nm. 
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Figure 12. Spectral temperature and fractional area determination for hot spot E. 



Table 1 . Hawaii hot spot temperature and fractional fill determinations. 
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Figure 13. AVIRIS measured spectrum, MODTRAN modeled solar reflected spectrum and MODTRAN 
modeled upwards transmittance spectrum for hot lava target C. 
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The AVIRIS imaging spectrometer measured the active volcanic vent region on the Big Island of Hawaii 
on April 26, 2000. Examination of the spectral image data set showed distinct expression of seven lava 




hot spots. These hot spots show enhanced radiance signatures in the spectral region beyond 1000 nm with 
respect to all other Earth environment components in the AVIRIS data set. The full spectral shape 
measured by AVIRIS provided a basis for temperature and fractional area estimation from these hot snot 



Anderson, G. P., J. Wang, and J. Chetwynd, “An Update and Recent Validations Against Airborne High 
Resolution Interferometer Measurements,” Summaries of the Fifth Annual JPL Airborne Earth Science 
Workshop, JPL Pub. 95-1, Vol. AVIRIS Workshop, Jet Propulsion Laboratory, Pasadena, California, 
pp. 5-8, January 23, 1995. 
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(AVIRIS), 1995,” Proc. Fifth Annual Airborne Earth Science Workshop, JPL Pub. 95-1, Jet Propulsion 
Laboratory, Pasadena, California, pp. 33-34, 1 995. 
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2919-2934,1993. 


ACKNOWLEDGEMENTS 


^ GO 

i.l 

a a 

^ (D ^ 

h 

(U o 
O O <D 

let: 

“ § ° 
3 

“ S 3 

c/3 d) c/) 

C/3 C 

O 4-» HH 

"irn'e 

£ 

(O 


goo 

i -2 

§ § J 

G> ^ C 

,2 S I 

4h O C/3 

T3 3 
Oh 

^ O 


^ o ^ 
o L- 


bps 


S O 


o 

00 



Measuring the Spectral Expression of Carbon Dioxide in 
the Solar Reflected Spectrum with AVIRIS 
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Figure 1. Plot of measured carbon dioxide concentration over time from Mauna Loa (Keeling and Whorf 
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EXPRESSION OF CARBON DIOXIDE IN THE AVIRIS SPECTRUM 
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Figure 3. A transmittance spectrum of the atmosphere with and without water vapor convolved to AVIRIS 
spectral resolution. 
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To explore the potential for AVIRIS to measure changes in carbon dioxide, an AVIRIS data set acquired 
over Pasadena, California was identified for analysis. This data set contains the lower elevation areas of 
Pasadena and the higher elevation areas of the San Gabriel Mountains. Figure 6 show a single-band 
image of this data set with Pasadena in the center and the San Gabriel Mountains to the left. The 20-by- 




spectrum. This simple analysis shows that calibrated AVIRIS spectra are sensitive to the differences in 
atmospheric carbon dioxide. The differences in this case were caused principally by the difference in 
path length through the atmosphere for the spectra analyzed. For this case, the AVIRIS spectral region of 
greatest sensitivity is indicated by the larger ratio values for the carbon dioxide band centered near 2010 
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Figure 7. Spectra extracted from Pasadena and Mt. Wilson regions of the AVIRIS data set acquired 
June 22, 1999. 
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Figure 1 1 . AVIRIS carbon dioxide image in atm cm. Lower carbon dioxide values are derived for the 
higher elevation regions to the left. Higher values of carbon dioxide are derived at lower elevation in the 
center and right of the image. 
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Exploring a Black Body Source as an Absolute Radiometric Calibration Standard and 
Comparison with a NIST Traced Lamp Standard 



O. 




AVIRIS radiometric calibration is based on an irradiance standard lamp and a calibrated reflectance 
panel. The lamp is a 1000-W quartz-halogen tungsten coiled-coil filament lamp purchased from 
Optronics Laboratories, Inc. The irradiance lamp is held a known distance (nominally 500 mm) from the 
panel in a fixture with baffles to reject scattered light. For radiometric calibration this fixture is 
positioned such that AVIRIS views the panel at 45 degrees inclination. The AVIRIS instrument response 
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Figure 5. Planck functions for a cavity black body at different temperatures. The data are scaled to a 
peak value of 100. The scaling factors are 0.00101 , 0.00763, and 0.258 for the 3000, 2000, and 1000-K 
spectra, respectively. 
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Figure 6. Results of fit between lamp irradiance and two-temperature Plank function calculations plus a 
fractional fill factor over the range from 1000 to 2500 nm. 



values delivered with the NIST traced lamp is from 2 to 6 percent in this spectral region. It is reassuring 



the NIST traced AVIRIS calibration fixture can predict the radiance of the 820-K cavity black body to 
this level of accuracy. 
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Figure 7. Planck function calculation for 820-K black body source In the AVIRIS calibration laboratory. 
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Following acquisition, the surface measurements were converted to surface reflectance using the 
interspersed reflectance standard measurements. The absolute reflectance of the reflectance standard and 
the corresponding bidirectional reflectance distribution function (BRDF) for the average 20-degree solar 



zenith angle were used for this reflectance calculation. Figure 1 shows the calculated average reflectance 
of the calibration target from the 890 surface measurements. Also shown is the standard deviation of the 
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from 4.1 to 3.8 precipitable mm. A value of 3.9 mm is estimated for the time of the AVIRJS over flight. 
This value, in conjunction with the surface reflectance and atmospheric optical depths, is the in situ 
parameter required to constrain the radiative transfer code to predict the radiance at the AVIRIS 
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Figure 4. Derived instantaneous water vapor abundance from the 940-nm channel of the solar 
radiometer data set. 
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Figure 8. AVIRIS average measured data for the calibration surface. 






predicted radiance for the calibration target on Rogers Dry Lake at 18:39:43 on June 5, 2000. This 
predicted radiance spectrum forms the b^is for validating the inflight calibration of AVIRIS for the 2000 
flight season. 
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to uncertainties. In addition, uncertainties in the in situ measurements and derived parameters will 
contribute to uncertainty in the MODTRAN predicted radiance. However, the primary result of this 
inflight calibration experiment is a baseline validation of AVIRIS inflight calibration at the 96 percent 
level for the 2000 flight season. 
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Figure 13. Comparison of the AVIRIS-measured radiance and the MODTRAN-predicted radiance from 
the inflight calibration experiment on June 5, 2000. 
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The standard deviation of these dark signal data provide an estimate of the instrument noise. Figure 15 
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corresponding values are given for 1999. 
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Annual JPL Airborne Geoscience Workshop, JPL Pub. 92-14, Vol. 1, Jet Propulsion Laboratory, 
Pasadena, California, pp. 1-4, 1992. 
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Oleoresin Chemistry and Spectral Reflectance in “Stressed” Lodgepole 
and White Bark Pine, Mammoth Mountain, California 



approach used with hyperspectral data has been principal components analysis (PCA). Rowan, et al. 
(2000), used PCA to identify absorption feature patterns obtained from vegetated areas with discrete 
bedrock geology or mineralization as the substrate. Many researchers highlight the need to advance our 
ability for hyperspectral imaging in vegetated areas as a near-term priority (Vincent, 1997 and Sabins, 
1999). 
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(Hickey, unpublished data). These observations support the premise that a discemable relationship exists 
between changes geochemical substrate and coincident changes in oleoresiri/monoterpene chemistry. 
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Each dataset was evaluated independently to facilitate feature selection and investigate any 
variance as a function of apparent plant health. An overview of the processing steps employed is shown in 
Figure 4. After the pertinent features were selected from each set, the covariance between the 
monoterpenes and spectral features was evaluated using ANOVA and regression techniques. 




features reported in Table 1 are derived from C-H bond variations prevalent in monoterpenes. It is 
important to note that other organics (e.g., sugars, waxes, etc.) also contain similar bonds. Further work is 
required to increase our understanding of the absorption feature source before a direct correlation between 
the absorption features and oleoresin/monoterpenes can be established. 
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higher level of complexity in the intermediate health classes. Identification and separation of the abiotic 
and biotic influences is, therefore important. Data processing and interpretation is continuing in this 
regard. 
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only to the health state but also the nature of the impact (i.e., biotic vs. abiotic). Included are the 
identification of a series of absorption features and monoterpene composition that discriminate between the 
two end-member health states (healthy vs. dying/dead). Additionally, covariance and regression analysis 
between the absorption features and monoterpenes shows that a statistically significant model exists that 
also separates the healthy from dying/dead end-members of the conifers. 
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The components are ranked in descending order of the Eigenvalues (image variances). The Eigenimages 
associated with large Eigenvalues contain most of the image variations while the Eigenimages associated with small 
Eigenvalues are usually noise-dominated. The majority of variations in the scene are generallv contained in the first 
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where 5, L and n are numbers of samples, lines and spectral bands, respectively in the image. The number of data 
points in the compressed image and the corresponding Eigenvectors are shown below. 
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Figure 2(a). Eigenvalues of the SCAR-B HSI data. 





Figure 3. Sample principal components of SCAR-B HSI data. 
Each component highlights various features in the scene, 

1^ PC: cloud, water, shadow; 2"^^ PC: vegetation; 

3^*^ PC: cloud structure and 4'^ PC: burned area. 
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Figure 4. Comparison of reconstructed and original data for sample features in the SCAR-B scene, 
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Figure 7. Comparison of NDVI derived from reconstructed and truth data over the Moffett Field scene. 
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the park such as north-facing slopes. Soils derived from rhyolite flows within the park have 
relatively low nutrient content; in these areas, the dominant forest type is lodgepole pine 
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5. MAPPING RESULTS 



altitude cover map in Figure 9. There is generally good agreement in mapping most of the 
as lodgepole pine cover. The low altitude data with its smaller pixel size shows a slightly 
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very little variation between the two data sets. This suggests that, despite the limitations of the 
combined radiative transfer and ground calibration method of atmospheric correction employed 



>■ ^ 


H M 


I ^ P 

2 S b 

Cl. • 

bO Jh >; 

Ml 


5 § • 

) O C 

3 *c/3 ^ 

3 ^ P 3 ’ 

J s 
> 

! I I ' 

i .3 

!-§| 

^ C'l S ^ 

rs «/3 S J 

(U a 

J ' 

^ w -o : 

»N p (D 

c3 t3 ^ 


^ i 

on ^ ' 


'O ^ 

(D bX)«: 
? § O 

p * s 
.2 gP-Sl 

.131 

^ ITi ^ ^ 
. On p kJ I 

^2: > S- •: 

c/3 ^ ^ 

r • ‘3 < 

H ^ § ii 

: o' M ^ ^ 

i life 

! I -S ^ ^ 

I ^ p 

! § § 

u ^ rt 

S 2 

S ic 
•HS |o 

Co lA 

.S ! 


" i 

, I 

: 

'O ® 

I ^ 


o^ 

o^ 


^ fe ' 

^ 2 s I 

(U o 
^ p <1^ 

^ <D tn, 

Tj C/) (Z 5 > 


246 



« &; 

s 

I ^ 

3 

I ^ 

cd Co 

0 J S 

" psi 'u ^ 

C- 

<D ^ ^ 

^ B ^ 
O I I 
■go -1 

§ o 

ffi a ^ 

fa Q .« 

<S-g 

' s2| 

r a> S 
' O > to 

io< I 

r r ftfl s 

; PQ .S -Si 

1 W ^ ^ 

jBre § 

; -g 13 ^ 

I U ^ 

' X ) F . 
r.s ,0 ^ , 

I (U ^ ^ 

x^ -h: 

1 a ; 2 ^ 

N '5 

i ^ s ^ 

o §: 

’<; g I 

I • 

■ p=i g 

^ U '* ^ 


ro O^ 
On ^ 
ON o 

csg, 

(U >. 

C 

5P s 

j2 Co 

"3 ^ 

. O 

_4. 

•> § 
d<^ 

*^0 Co 
ro 

cd 

o CO g 

•Jl 3 rN c/) , 
cd ^ p , 

X) o 
P Oh' 
’ Ph CO 
cd 


g^'g 

P ci. 

§:<< it X 

sm 


^ §-^ 
fai w s 

T3 -d El 

§ g I 

O'® * 


- is >5 

« u "ts 

X W) Oh S“ 

2?.S w i; 


^ Q i ^S 
-■'' '-' '3 -a 

(D O [C 

+-* faj 1 ^ 


'oh S 

PhX 
C5 •'^ 

S ^ 

3 X tjo 

■S --s .s 

EC ^ bo 

fa ^ 
3 M t 

QQ S h 

CM 3 « 

fl CT'cfc- 

.2 ^ 

I3 ' 

Ph 'S ^ J 

S cd ^ 

O (U ’ 

^ ^ V 

< : 
'Td cd ^ 
.2 X s • 

S o ■ 


I ¥ 

<H 

<D 

>• 

•p ^ 

• cd 

NO .ti C^ 

^ A ^ 

t^ O' 


>;5 ^ 

S X t 
o 3 s 

Os U O 

.c g 


< d -b 

r)A ^ O 'X 

O cn w 

r bfl S S 

^ -s I ;« 

pd g ^ 

H rA> 


I 

SS kV 
. X ^ 

21^ 

X) 5 Cj 

g Q g 

. M 3 
CO 2 « 

I rP rt 

<N <DJ} g P 
- iS P rrf 
00 X is 2 
(N Cd o 

c^ C/5 fa 

S <r' >5* 
•2^..S 8 
^ SP g 

o Z g £ 

gif's <2 


O « 

2 5 
o 8 § 

O C/) ^ 

is S I 
S. §1 

'M Cd 
b-J cj CJ 

pd ^ 

S fa 

ci, 9^ <3 

S^" HH BO 

^ S 

b 2 P 

^ c/3 (D 

t § g 

§ c^ g 


^ "g - 

'O S 2 

g . ss 

I Z -S’ 

^Sl' 

^ O ' 


<u 

> 

S ‘p 

fg g 

Id CO 

o S 

CJ P 


® 2 l.s 

MSI 

^2 5!^ 

« § 2 : -I 

o ^ .H s 

.. 2 w 


> Bo^b 

I S §Q 


Ktsi 

o P S 

ON P O 

^ -S GO 
r 4 P 3 o 

rin <U (L> O 

^ ^ ffi § 


& O 

2 

S "2 

5 X 

5- CO 


o SP 

« .S 

c/3 ^ 

’go ^ 


P 05 
(D S 

•8 feb 

O <L> 


r^ <D 

O >H 
0 H-S 3 
W .2 S 

’^^•1 I 

O H ‘i:* 
> P ^ 


O rf P y fS 
2 ; •§ 2 ffl « 
o'* 

(U d 


„ Cj 

; fa 

S g 

hT «3 3 
^ dx 

rv S ^ 

' a-s 


3 ca Bo 

-s 

^ 2 § 

' Pi .2 

T S ^ 
_ 'o 

h 't:? 

■6 § 1 

S 02 

•S I e 

CB *S X 

Oh 2 S 
“ (L) X 
^ BO 3 

Q C /3 

fa W)cO 
i <=*^ faS 

1 ^ g *> 

' ^ s c 

tS 


iT' 

(D > 

4-> ^ 

O 

XJ o 

O c/3 

tl 

2 c/3 

S g 

O S 
cd 

00 ^ 
ON .2 

^ o 
Os & 

C2 o 

C /3 


. '^’SP 

^ p 

_• r^ bD 


13 

K |V 
O ^ 
•b M 

s M 

y U 

^ r 
' <«>fa 

' .S Pi 


r- 

'd- 

<N 



Martin, M.E., Newman, S.D., Aber, J.D., and Congalton, R.G., 1998, “Determining forest 
species composition using high spectral resolution remote sensing data,” Remote Sensing of 




O Z 
fti a 


PI : 

2 r <u 
S O <u 

2 tT ,’Sh o< 

§ ^ 3 ^ ' 

00 -S MO* 

r 1 g ^ 

-1 c ^ 

^ I .s I : 
5 1 1 ' 

^ I 

• o ,e 
^ U « ^ i 

£, H ,5 ! 

: S CO 2 

; D CO Co ^ 

i r-Sa- 

I 3 ' ^ ' 

’ Q fcl i2 
; . 2 « ’ 
(/3 ( 

i t: ^ o 
I 2 o S . 

O tJ 

Pi 


wCo § 
§ 

2 '“ >> 
> -S 

iti 

C<3 ^ 

«f ^ W 
(DON 

2 S 2 

s 

c . 
> « ^ 

"S' s ^ 

5 


2 I. « 00 

M C 

^ .2 ^ ^ 

h 4 iS 'g ^ 
“ (U o 
00 OD § 

•S > •< s 


1316. 



248 



Table 1. AVBR.IS Data Collection Parameters for High and Low Altitude Acquisitions. 
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spruce/subalpine fir, SF) system as the best match. 
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Forest Cover Map from High Altitude 1996 AVIRIS Data 
Old Faithful Area, Yellowstone National Park 
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reCambnan bedrock m the NGM area consists of limestones, dolomites, sandstones and their contact metamorphic 
equivalents, however, published geologic maps do not distinguish between the different lithologies. Mesozoic 

mapped primanly as granitic-composition and some age-dates are available (Albers and Stewart 
iyll). Mesozoic units mannpH in : - i i ’ 
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pervasive argillic alteration that is characterized by the formation of alunite and sericite (muscovite); 




. spatially restricted, phyllic wall rock alteration which is characterized by the presence of illite. This alteration 
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The Mineral Mappi ng Process: AIG plans to apply the following well-extablished approach to mine ral map p in g 
using hyperspectral data to the Hyperion data. These methods have been tested extensively using AVIRIS and other 
airborne hyperspectral datasets (Boardman and Kruse, 1994; Boardman et al. 1995; Kruse and Huntington 1996; 



Figure 2: The AIG “Hourglass” method for analysis of hyperspectral data. The underlying image is the 
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contributions from carbon dioxide, ozone, and other gases (Gao and Goetz, 1990). Strong atmospheric 
water absorption bands make the atmosphere opaque in many regions (for example the 1.4 and 1.9 jam 
regions) and only small atmospheric windows are available for terrestrial remote sensing. 
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identified using a spectral library (right) and Matched Filter mapping results (left). Brighter 
pixels in the Matched Filter results represent higher abundances. 
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the utility of Hyperion for geologic applications. 
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HYPERSPECTRAL IMAGE CLASSIFICATION USING A SELF-ORGANIZING MAP 
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by synapses or connecting links, each of them characterized by a weight. Specifically, a signal xj at the input 
of synapse j connected to neuron k is multiplied by the synaptic weight w^j. 



Figure 2. SOM neural network topology including weight matrix W and self-feedback and lateral feedback 
connections. 
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d) Learning: adjust the synaptic weight vectors of all neurons, using the update formula (4), where ri(t) is a 



5. RESULTS AND DISCUSSION 
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ground truth classes, associating exactly one neuron with one class. Usually this fact is not possible if we 
have a larger number of output neurons that ground truth classes. In this way, one correct classification 
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Table 1. Resulting confusion matrix for the Indian Pines dataset considering 100 iterations, 16 neurons in the 
output layer and Gq = 2. Non-zero cells in the matrix are highlighted. 


In our second experiment, we consider a smaller number of neurons in the output layer (5), 80 iterations and 
Co = 4 . Figure 4 shows the resulting classification and the associated confusion matrix. We can appreciate an 

improvement in the topology of the confusion matrix (horizontal “lines” are more contiguous). In this case, 
the overall accuracy is 60%, without considering the unclassified pixels. 



.2 ^ 


S) 




•a ^ 












(U 

73 

c 

o 

o 

§ 

o 

(U 


T3 

d> 


I 

C/5 


(L) 


° I 
B I 

c3 on 

Cii ^ 
C cd 


o 

(L> 

CU 


O o 


<D i/J 

IS ^ 

<l> 

t.s 

O ^ 
cd bX) 


cd (D 

S 

3 2 

■o 

O 1) 

o 

M 

c« 


C w o 

^ t/) *r 2 


I 

cd 


o .22 


13 

O 


Ph 

(L> 

6 

a 


^ §: 

H S' 


<N 

r-- 

<N 


It is more efficient, in terms of accuracy, to increase the number of neurons and decrease the number of 
iterations. As we increase the number of neurons, the overall accuracy increases, but the topology of the 
classes is poor. 
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Figure 5. Resulting classification and greyscale representation of the confusion matrix for the Indian Pines 
dataset considering 20 iterations, 16 neurons in the output layer and Gq =16. The time used for this 



Table 2. Other experiments performed over the Indian Pines dataset. 
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distributions of alunite and pyrophyllite. 
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magenta pixels show this 0.62pm feature, and hence likely contain cyanobacteria mats in some proportion. 1 992 
AVIRIS data is similar to that data expected from such sensors as Hyperion. Hence these results suggest that simple 
vegetation zonation patterns and communities of thermophilic organisms indicative of hydrothermal discharge may 
be identifiable from space-based hyperspectral sensors. 
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three datasets, although AVIRIS (17-1 8m) has a poorer spatial resolution than HyMap (3-5 m). Finally, alunite 
mineralization is detected in both AVIRIS and HyMap, although some deposits were ultimately missed by AVIRIS 
due to its poorer spatial resolution. 



SNR is obviously a controlling factor for success in mapping, however the “poor” SNR of 1992 AVIRIS 
data did not totally hinder mapping efforts. In fact, the mapping success achieved in 1992 data suggests optimism 
for space-based hyperspectral mapping efforts. Additionally, spatial resolution is clearly important for biological 
mapping, and can be critical for some geological applications. 
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Figure 2. Distributions of mineralization at Rhyolite Hot Springs in 
the west-central caldera. Glassification done in ENVI with Matched 
Filtering. 
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Figure 5. Matched Filter results for AVIRIS and HyMap imagery. Mapped alteration in 
yeliow is alunite. Although the AVIRIS data does reveal alunite, the full distribution is 
only mapped by HyMap. This is likely due to differences in spatial resolution {17-1 8m 
for AVIRIS vs. 3-5m for HyMap). 








DECISION TREE CLASSIFICATION OF MULTIDATE AVIRIS DATA FOR MAPPING 
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ENVI is fast and the method is intuitive. Endmembers are readily identified, but actually 
determining the mineralogy within each pixel is not easy even using ENVI’s built-in expert svstem 
the Spectral Analyst. 

The results of applying ENVI to the CF AVIRIS imagery are shown in Figure 3 The 
primary endmember spectra plots are shown in Figure 3. 


■o = 

0 S^ 
.bz 0 
0 ^ 
0 0 O 

■°#Q. 
0 . 
S 52 o 


o 

0 0c 

•“ ^ .9 

•= 0 .0 
-C Q. *43 

I— 0 CD 


6 g»£ «3 

: a g 

^ I- E (0 

1 CD IS 

3 ds •“ 45 

L ® S ^ 
2 -c 

3) 0 CD 

is = 
: i-§ 5 

'■U° 

12 0 0 
> 0 Q. 0 

• T 3 ^ . 

* 0 5 ' 

t jC £ T) 

: c 

c ^ o 
' 0 c^ 

i 0 p 0 
-.Sow 

W 5 o 
Q. C 
;E CD Q. 
0 -5 C 
C ^ 0 


292 





m 

ON 

(N 


absorption feature identification. 
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Discussion 

The results presented here are meant to highlight the primary techniques for processing 


S « 

■D ^ 

<0 y ^ 

CD — 

> -J= > 

2 ^ Z 

til LU 

>'E • 

T3 — CD 

0 > -O 

“ 7 ' — 

2 

c" m ^ 

E 2 0 

CD 0 Ql 

0 C 2 
0 0 0 

3 O) a. . 

c x: 

-C . O) P 

1 gi g, 

^ Sit 

■D > CD a> 

2 <C ^"D 
0 0 .E I- 
S ^ CO 0 
S ^ 0 -D 
0^0 1 - 

CD o C^ 

> > c: 0 

^ O) 

0 C O ^ 

m O h- 
0 0 ^ 0 ) 
jc :e 
/) <1^ a ^ 

CL CD ^ 

E. o O 

0 0 0 -^ 
CD O .f_i E 

2 CD s ® 
E ^ P F 

3 S “ “ 
g - = o 

§■8 g >> 

CD E £ > 

C^ ^ CD 
O £ 



Determine 
Pixel Content 


Figure 7 Hyperspectrai data processing flow showing where the respective software 
process is introduced to the data. 
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SOM are strictly a classification method that does nothing regarding absolute 
interpretation of pixel content. The result of using SOM is that larger aggregates of minerals, i.e., 
rocks, are discernable in HSI data. SOM represent another classification paradigm offering 
significant benefits for unsupervised image classification. 
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The AVIRIS data system is described by Green, et al. (1998). One AVIRIS flight line was collected on July 
20, 1999 (14:00 GMT) in the New Bem-Kinston stretch of the Neuse River, immediately upstream of the Neuse 
Estuary (Fig. 1). Aircraft motion geometric distortions are minimal. Minor haze and clouds (highlighted in images at 
1882 nm) are considered negligible for this analysis. The single scene analyzed here 

(f990720t01p02_r02_sc06.vl.img) is dominated by agriculture and wetland forest flood plain, and is almost entirely 
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2.2 Endmember LCLU Linkages 

A recently compiled (U.S. EPA), remote-sensing based Neuse River Basin LCLU map was used to identify 
the land cover class of pixels identified as AVIRIS spectral image endmembers. The map was produced from 1998- 
99 SPOT 4 XS and Landsat 7 ETM* scenes, plus CIS input data (roads and wetlands). It has a pixel size of 15 m, 
with a 0.4 ha minimum mapping unit (MMU) (0. 1 ha MMU in riparian zones). Overall Level I accuracy as verified 




The Virtual Field Reference Database is available online at www.epa.gov/nerlesd 1/lcb/nrbA^FRDB/ 
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Figure 5. AVIRIS scene f990720t01p02_r02_sc06.vl.img and locations of image endmember pixels. Above: 
AVIRJS image (visible bands) with endmember pixels plotted. Below: endmember pixels only (image removed). 
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Inspection of spectra at 650 nm shows the endmembers to break into three groups (Fig. 4): 
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SPATIAL/SPECTRAL IDENTIFICATION OF ENDMEMBERS FROM AVIRIS DATA 
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The extension of classic morphology to the hyperspectral domain allows us to integrate spectral and spatial 
information in the analysis process. Tn Sertinn ^ Crimea _ , • 1 , . ^ 
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The principles of mathematical morphology have also been extended to the grayscale image case (Sternberg 
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of any family ofN-dimensional vectors. 


The definition of a vector ordering relation can be done as follows. Let Xi, X 2 , x„ be the hyperspectral 
image pixels within a filtering window that represents the structuring element of a morphological operation. 
We can define a measure of the dissimilarity between two of those pixels, Xa and Xb by: 
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Output of the operator=dist(V, E) 
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3.3.3 Redundant endmember thinning 
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Figure 4. Intermediate results of processing the Salinas A image with AMEE method, a) Probability image 
after applying morphological operators to the full hyperspectral dataset, b) Resulting image after 
automatically thresholding a with Otsu method, c) Resulting image after spatiaEspectral growing of the seeds 
obtained in b. 
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We have presented a new autonomous method to find endmembers in hyperspectral remote sensing data 
based on mathematical morphology concepts. Although no previous dimensionality reduction or data thinning 
is performed, it is capable of accurate and rapid determination of endmember spectra from high dimensional 
hyperspectral images. Some preliminary results obtained after applying this algorithm to AVIRIS sensor data 
show that it is successful in the task of integrating spatial and spectral information in the analysis process. 
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Szeredi, T., Staenz, K. & Neville, R., “Automated Endmembers Selection: Part I Theory,” Remote Sensine of 
Environment, 1999. 

^ _ 


S P S:; 
^ P 

U 



318 




Figure 5. A subscene of the Salinas 98 dataset, showing known ground truth as polygons containing solid 
colors bounded by black lines and surrounded by thumbnail figures of the ground truth taken at the time of the 
data acquisition. The outlined box in red shows the subscene boundaries of the scene called Salinas A. 
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DETERMINING ATMOSPHERIC COLUMN WATER VAPOUR IN THE 0.4 - 2.5^im 

SPECTRAL REGION 
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L{W,X)= taMW" 



The following procedure is then implemented to estimate the atmospheric column water vapor W, 
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synthetic signal that is used in step 1 , and Table 1 has its limits changed to those given in Table 2. 



modelling program, such as MODTRAN4.0 (Berk et al, 1999; Adler-Golden et al, 1999) to produce a series of 
synthetic at-sensor radiances that differ only in the amount of atmospheric water vapor. For the purpose of the 
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where, it is assumed that the solar zenith angle //„ and the sensor zenith angle n are known. Further, U is assumed 
that the spectral solar irradiance £■„(>!) at wavelength A is the same in the atmospheric modelling program of choice 
and the actual spectral solar irradiance observed on the day of acquisition. (A) is a synthetically generated 
estimate of the path radiance, and is given as, 



<N 

m 


is preferable to weight consecutive absorption features (ie the measurement channels). This can be thought of as a 
transmittance slope and is analogous with the reflectance slope that is used in other methods. The weighted mean 
radiance is given as, 



To estimate the atmospheric column water vapor W for a given wavelength >i,.,.the change in scaled 

radiance as a function of water vapor needs to be established. Writing the scaled at-sensor radiance of equation 
( 1 0) in a power series form yields the following, 
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evaluate the accuracy of the procedure have shown an accuracy of around 1-2%, in reality though this is expected to 
be higher due to incorrect modelling and other factors such as aerosol contamination and random instrument effects. 
It is expected that this procedure can be used in a Look Up Table (LUT) procedure. A comprehensive sensitivity 
analysis will also be carried out in the near future. 
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Atm o spheric Column Water Vapor (cm) 

Figure 1: The radiance at-sensor as a function of atmospheric column water vapor and for a given surface spectral 
reflectance. Four different wavelengths are shown, they demonstrate that the radiance can be represented by an 
approximating polynomial. The order of which is dependent on the region under consideration. As the water vapor 
absorption feature (940nm) is approached the order of the polynomial is seen to increase. 
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These studies were conducted only with funds and/or resources committed hy the participating 
bureaus. Studies were selected on the basis that baseline information pertinent to each of the 
studies proposed was already well developed and readily available. NASA contributed AVIRIS 
missions for each of the selected sites, and provided technical assistance both directly and 
through partnering organizations. 
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charactenzahon of hyperspectral data, 3) exploitation and expectations of the data, and 4) data 
fusion techniques, degrees of difficulty encountered, and known successes and failures in the use 
of the data and/or technology. Unfortunately funding limitations did not permit the 
implementation of this last element of the technology transfer project. 



A siunmary of each of the four study sites is provided below including the personnel involved, 
objectives, and brief background information. Accomplishments, benefits to DOI land managers 
and study site partners, and lessons learned are discussed. 
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species and vegetation associations. However, multidirectional geometric anomalies of up to 
2 to 3 pixel widths, which were formd in both the low and high altitude data are inhibiting 
use of the data for mapping purposes. Both scenes were formd to contain non-systematic 
geometric characteristics that cannot be eliminated via traditional registration techniques. 



Study site 2 - Estimating the effect of invasive woody species on grasslands, Great Plains Basin. 
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reservoir watershed management alternatives to limit the amount of Mercury influent to the 
reservoir and to venfy the validity of making direct remote measurements of water quality in 
the reservoir, m terms of total suspended solids and chlorophyll content. 
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A significant outcome of this study was BOR’s decision to acquire a full-range field 
spectrometer after seeing the demonstrated potential of imaging spectroscopy for assisting 



The primary overall result of the NASA/DOI hyperspectral technical transfer project was the 
successful demonstration of AVIRIS hyperspectral imaging capabilities to land managing 
bureaus within the Department of the Interior for: 

( 1 ) Detection/mapping of multiple types of invasive plants 



participating DOI bureaus, NASA, the USGS, and academic partners. 



Unexpected outcomes and benefits 
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connections between individual scientists and bureau disciplines were enhanced. Participants 
from all projects made ongoing commitments to the technology ranging from obtaining 
hyperspectral imaging software to procuring entire ground and aerial data acquisition systems. 
In two instances, study participants are continuing their work through partnering efforts, and by 
providing their own resources and support. 
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MODO: AN INTERFACE TO MODTRAN FOR THE SIMULATION OF 
IMAGING SPECTROMETRY AT-SENSOR SIGNALS 
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The core interface of the procedure is the tapeS editor window. It allows to set most of the input parameters using 
pull-down menus instead of manually editing the rigidly formatted ASCII file. Logics within the tapeS are considered, 
such that if, e.g., the transmittance mode has been selected it is not possible to set the irradiance source options. Sub- 
interfaces will pop up for supported special functions such as the import of user defined atmospheres, the selection of 
the surface reflectance, or the definition of the four standard aerosol layers. The interface is grouped in the same way 



MODO 

input preparation 
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Figure 2: Example of the MODTRAN interface ‘MODO’; translator and organizer for the ASCII input file ‘tapeS’. 
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Figure 3: Interface for the selection of spectra from the spectral albedo file for the MODTR AN background 

reflectance option. 



The typical workflow using the MODO utility depends on the task to be performed. In imaging spectroscopy, the nor 
mal case starts with known surface reflectance spectra which need to be transposed to at-sensor radiance values. For 
the creation of spectral databases or Look-Up-Tables for later inversion, standard setting for reflectance and discrete 

values for parameter variation are taken as basic input. An overview of a typical data simulation workflow is given in 
Figure 4. ^ 




n 


Figure 4: MODO typical workflow structure as used for the simulation of imaging spectrometer data. 





The standard unit in cm~^ is given as the original MODTRAN wavenumber reference which may be related closely 
to the energy levels of the simulated photons. In imaging spectrometry and spectroscopy of the visible/near infrared 
part of the spectrum, the usually used wavelength reference are microns or nanometers. As the resolution of current 
imaging spectrometers is in the range of 5 to 20 nanometers, it has been decided to select the wavelength in nanome- 
ters as generic reference for data simulation within MODO. 





convolved data values in any case. 
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MODTRAN 4 Radiance 



Figure 7: Total at-sensor radiance over vegetation convolved to specifications of the APEX instrument and to 

reference instruments (AVIRIS, HYMAP, and DAIS 7915). 



t)X) 

•S - 

c 3 Cd 

E o 

if 


S X 

o W 

o 


cd 


E 

^ S 

<u ^ 

Cl 

^ P 
(U <D 

-C ^ 


5=! 

^ O 

(D O 
v-i cn 
3 ^ 

8 


Uh 

c ^ 

•!=! C 

c .2 

<D 3 

-E 

’51) 'o 

^ JH 

cx .E 

§ s 

cd 0) 


O 

cx 


X 
<D 

^ 2 

c "O 

E (U 
o J!; 
•n cd 

a 


° I 

o 


o 

C/D 

Si 






o ^ 

a 2 

So <u 

E E 


fX 


o 

(D 


T3 

C 

cd 


§■.1 

S O' 

4-1 (L> 

cd c 
c/) p 

*3 c5 
_§) o^ 

C/D ON 

ON 


Id 

2 

(D 

j:: 


<D 

3 
3 

1 o 

S ^ 

p w 
^ 00 

(/> C/D 

3 3 

e I 

2 “ 

(D J2 

1.0 

o S 

Cd 

o 


x: 

cx 

o 


o 
c 

CZ) 

8 

TJ 

i-C (]^ 
O 3 

Cd o 

2 E 
a . 

9 - T 3 

3 (U 

w 

T3 cd 

B 3 

S 2 

SS « 

E ^ 
CX X) 

2 O 

i 

S S 3 

s s 

E 


o 

X 


ifl 

Z 

o 

hN 

C/) 

nJ 

u 

z 

o 

u 


uo 


(L) 

•E 
c 


E ^ 


ti .S 3 

|o 

^8 


W) 

c 

'§* 


.'e 

3 

o 

T3 

CD 

C; 

B 

C/D 

CD 

CD 

2 

s 

C/D 

CD 

E 

2 

00 

X 

3 

c^ 

CD 

Td 

<D 

E 

2 

.E 

•C 

O 

C/D 

c 

o 

CD 

X 

cd 

O 

o 

2 

E 

C/) 

•D 

T3 

C 

(D 

CD 

X 

O 

X 

cd 

c 

2 

Vh 

1 

X 

?> 

O 

X 


J 

CD 

cd 

Oh 

O 

Wh 

t2 

r-i 

s 

2 

X 

2 

C/5 

C 

o 

ling 

< 

S 


<u 

o 


Id 

b 

t: 

2 

c/5 

2 

o 


S 

2 

u. 

3 

6 

3 

CD 

’"3 

2 

E 

2 

Q 

*E 


(U 

o 

2 

(£ 


&* 'oi 2 

E 'H & 

•-oi r-< tJ 




CZJ 

H 

Z 

U 

U 

O 

Q 

X 

o 

z 

X 

u 

< 


0 
00 

c/5 

1 

■3 

Id 

1- 

eg' g 

(L) 


C 

(D 

I. 

O 


. (L) 

X: T3 
(D 

^ X! 

3 ts 
N E 
^ c2 

o ^ 

^ 5 

.3 ex 
E C 

c ■§ 

L5 3 

<D ^ 

- ^ 

IS 

s > 

■c £ 

2 a 

2 o 

S > 

I— 5 D 

g>=5 

•i 

C t -4 
(D O 
00 ^ 
^ Oh 

0 T 3 

1 § 
o I 

X E 
^ ex 
>» 2 
X g 

T3 O 

0 CD 

ex 3 
ex X 

3 O 

C/D X 

c 

CD ». 

^ 3 

C/D 

J§ .E 

J= ‘o' 

■g « 

1 § 

^ J5 

</5 ex 

£ - 
H 00 


C /1 

U 

U 

Z 

U 

u 

U 



o 

»T) 

CO 


4 

X 



Characteristics of AVIRIS Band Measurements in Desert Agroecosystems in the Area of 

Blythe, California: I. Studies of Cotton Spectra 
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MATERIALS AND METHODS 
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nrst-order soil background vanations (Huete, 1988), transformed soil adjusted (TSAVI) (Baret et al 
1989), ATSAVI (Baret and Guyot, 1991), and water band (WBI). They were calculated as follows: ’ 
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gave positive correlation coefficients (r^ = 0.80- 0.87). Crop identification was explored using Spectral 
Angle Mapper classification (SAM), spectral feature fitting (SFF) and spectral binary encoding (SPE). 
Excellent agreement as found between the predicted and the actual crop type (i.e., the correlation is 
between 85-90% match). Additional analysis of the healthy cotton showed six local maxima. 
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Fig. 5. Healthy cotton peaks that were collected by the ASD FieldSpec®. 



r'^2=0.612254 SE=0. 176314 F=19.1338 



Fig. 7. Data extracted from AVIRIS image on relatively stressed cotton. 





r'^2=0.77082 F=355.005 




Fig. 9. Peak relationship between soil moisture and reflectance in healthy cotton. 




Pk=Lorentz Amp 1 Peaks Bg=Linear 
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Overall, the measured spectral reflectance signatures showed variation associated with different crops at 
different stages of development. This supports the continued research into use of reflectance spectra in 
the AVIRIS spectral range for agricultural applications. Future efforts will be directed to extend this 
research to other crop types and agricultural regimes. 
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Preliminary Results of Hydrothermal Alteration Assemblage 
Classification in Aurora and Bodie Mining Districts, Nevada and 
California, with Airborne Hyperspectral Data 
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Both high and low altitude data sets were ATREM corrected, but only the high 
altitude set was EFFORT polished; the procedure has created significant artifacts in the 
low altitude set. For the preliminary analysis, both sets were subsetted to SWIR range and 
individual scenes mosaicked. A Minimum Noise Fraction (MNF) transformation was 




measurements for classification. 
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supergene, clay in the tuffaceous sediments ringing the Martinez Hill. In Bodie, it is 
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Figure 1 - Area Location 



Figure 2 - AVIRIS and LAB Spectra 
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Aurora, NV - SWIR Spectra of XRD Samples 
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Figure 3 - XRD Spectra and ASD Field Spectra on XRD samples 
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Spectral Dimensionality and Scale of Urban Radiance 

Christopher Small 
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scale dependent optical properties of existing urban settlements may influence future design decisions. 
Charactenzation of urban spectral reflectance serves both of these objectives. 


%-* 

cd 

N 

2 

T3 

o 

o 

*-♦-» 

(Z> 

4-> 

o 

2 

C/) 

CO 

cd 

'C 

<D 

oi 

'C 

(D 

-S 

73 

(D 

O 

.2 

[-C» 

-4-^ 

o 

(D 

X) 

*C 

cd 

o 

13 

'bb 

o 

c/i 

J3 

4h 

-3 


(D 

o 

O 

o 

‘co 

73 


(D 9 f? 

e hi ^ 
o <D ^ 

~ 3 s a 

2 Cl 
O ^ 

|g Ss 

y Si s 

^ § Si 

s §3 g ^ 
9 o p S 
2 g € ■?< 

S o. ® 

g ^ 3 ^ 

2 o -.3 o 

0 'J- (D o 

^ .g -c 

=: S 2 

^ i -c .2, 

2 o '*5 2 

1 ■s ° <=^ 

^lll, 
«| 2 I ‘ 

^ 2 -S -S 

J 53 ^ -3 

^ S § 

^ o ^ 

' ^ ■§) ° 

3 9 p (u 
3 .2 o 3 

r c« (u 

: (D »H O 

^ T3 O ’ZJ 
i 73 6 

<D 9 

r Si ^ S 

i cd c« « 

; § 

» ^ &•'§ 
^.2 2i « 

I ^-2 6 

! 4h ^-t 

> PhH O 


iM 

^ ^ 

1 2 

ifl'2 

TTJ ^2 
^ C/i +-» 

<D (U <L> 

T3 O 3 
"5 'U o 

g « 2 

2 3 tj 

<L) Cl. 
O c« 

S.O 2? 

^ c/i C 

W) 2'S) 

•S -s § 

•3 

g 

-3 “ ii 
(D ^ 
^ 4h cd 
O ^ > 

§1 S 

<D 

3 “ 

^ § s 

fl) C cd 
fX rt 


i> g -s 
e5 -2 & 

^ 2 S 

•3 ° ^ 

gw© 
2 « 2 
© ^ § 

2 §‘^ 
O (D c/) 
_o 73 
fH 4-J cd 

q 2 

^ 

^ 7? cd 

§ |s 

£ 5 60 
^ 2 S 

2-- i 

.S.2 I 



VO 

m 


gives an estimate of 





a 


o 

o 


U 


Figure 1. AVIRIS false color composite of Milpitas California acquired in 1997. The Red, Green 
and Blue layers correspond to the 2.2 J {x/n, 0.8 J jiwi and 0.66 fiw wavelength bands respectively. 
The image is 6x8 fan and each pixel is approximately 20x20 m. The spectral dimensionality of the 
mosaic of urban, subriiban and undeveloped areas varies at scales larger and smaller than the 
Ground Instantaneous Field Of View (GIFOV) of the AVIRIS sensor. 
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rapidly diminishing low-order eigenvalues and the gradually diminishing tail of higher order eigenvalues. 
The successive reduction in the number of eigenvalues greater than unity is a direct consequence of the 
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Figure 6. Spectral resampling and loss of dimensionality. Resatnpling AVIRIS spectra hy factors of !/n 
for tt= 2. 4, 8. and 16 rapidly reduces spectral dimensiomlity. Resampling by a factor of 1/4 results in 
some loss of significant dimensions as indicated by the leftward shift of the breakpoint in the eigenvalue 
distribution. Resampling hy factors of 1/8 and 1/16 results in greater losses of spectral endmcmber 
resolution with spatially coherent eigenimages for all but the highest one or two dimensions. 



with the highest spatial autocorrelations are associated with the hyperspectral equivalents of the spectral 



endmembers that characterize broadband imagery of Milpitas and 
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The main objective of this study is to assess ISODATA classification for mapping land cover in Northeast 
YNP using select bands of low-altitude AVIRIS data. A secondary, more long-term objective is to assess the 
potential for improving ISODATA-based classification of land cover through use of principal components analysis 
and minimum noise fraction (MNF) techniques. This paper will primarily report on work regarding the primary 
research objective. 
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Figure 1. Location of study area within Yellowstone National Park. The graphic shows a hilishaded 90-meter 
DEM from the USGS with the AVIRIS scene boundary outlined in red. 

The number of bands clearly affects the amount of time needed to perform ISODATA classifications. 

Based on first-hand experience, the amount of time needed to run an ISODATA classification on a whole AVIRIS 
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Table 3. Bands selected for ISODATA classification. 
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Figure 2. AVIRIS-based final classification with 33 noiap categories. 



Table 4. Summary area for map categories in final classification. 



'BBW denotes grassland community with bluebunch wheatgrass 































Figure 3. Visual comparison of final classification to scanned 1:8,000 color infrared aerial photographs. 
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It is quite possible that there are other comparable or even better ways to compute wall-to-wall land cover 
maps from AVIRIS data. For example, some preliminary work by the author indicates that use of the 15 most signal- 



Adams, M. L., W. D. Philpot, and W. A. Norvell, 1999, “Yellowness Index: An Application of Spectral 
Second Derivatives to Estimate Chlorosis of Leaves in Stressed Vegetation,” Int. J. Remote Sens., 20(18):3663- 
3675. 
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3. Carbon estimates with multispectral coarse resolution remote sensing data 
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Leaf-level optical indices and ratios centered at the water absorption bands at 940nm and 1200nm show 
good correlations with ground-based measures of leaf water thickness (Ustin et aL, 1998; Gao and Goetz, 1995; 
Serano et aL, 2000) and therefore such indices can be used to map canopy water content. In contrast to NDVI, LAI 
estimated from canopy water content, remains linear to LAI > 10 (Roberts et aL, 1998). Estimates of canopy water 





•S < C 
° S ^ <1^ 


O, O 

H ^ 


»> in »-' Ch 

"3 5 O 


c-s^ s 

<D wj .ti 
« h 

^s i 


I bo D ^ 
■U > 2 

' 'S ^ ^ 

I r-i ^ 

2 §1 
D dn 
(/j O 


2^ 2 Tf 

3^0 


-i-t V—' (U , 

^ <D rH C ^ 

d CL ^ Cd 


■3 S 2 


^ « HH rrt 

I* g 2 “ 


g 

S ^ fe I a I 
|e^|| & 

<D t-H c^ a c/3 

T3 < ^ te ^ 


rvT' 2i 'u ^ 2 ^ 

B 3 T3 " g g 

> .2P 3 g § o 3 

^ § “ ^<N -S 

^§ 3 0 0 S ^ 

^ -a g* 3, S ^ ^ 

I 'S ^„.S^ ^1 

I III I ^ S 

g > ^ (D d O 

B <1^ "2 S 5 *5 (D 
i ^ g I ^ 3 S 

^ o g o cH r^ S3 
^ o g > -S tj ^ 

•S H Cr ^ ^ 

O bO (D S? 

d • ^ 9^ a ^ 2 


c/3 ^ S jn 

c w § o 


_0 ^ rrt 


“ 3r>:3 g S 

g.1>l 0 3 


'O ^ 

o 2 CO 
2 ^ (D 

§i§ 


3 «> S 

.3 3 3 


S3 C g 

111 S c 
o o 

<N 3 


“ £ s I 

b 4 S 
c/3 CJ ^ a 

s|^|S 

lllii 




%) 0OUB}O9|iay AdouBQ 


a w d 

- § ^ 
O O 1— 5 

^ ^ 2 

Cl 9 a 

w i5 43 


^ o2 

(/) '' c- 

•> 


E 

<D 

> 

2 

1 


(D 

CO 

>-l 

o 


d 

od 

c 

S) 

d 

<D 

'♦><» 

;§ 

CO 

SI 
•1— • 

CL 

4-J 

g 

E 

td 

i-i 

(D 

Id 

d 

Q 

D) 


CO 


a 

d 

}-H 


C 

0 

bO 

•n 

C! 

u 

o 

o 

C/3 

C 

o 

(D 

a 

-♦-> 


u 

a 

•a 

I' 

a 

3? 

tc 

CL 

CO 

<D 

I 

> 


a 

> 

T3 

d 

W) 

r2 

& 

3 


fc^o 

(-1 _9 

O H 


■2 X <• 
to D o 


3 3 T3 _ 
3 ^ 


UV 2 


- 3 .2 

~ "S 

« 3 .S 

■Jr* ^ I— < 
2 1^ 
3 S > 


^ Y d (D ^ 

X) ^ 2 ^ H-l 

? (D a <C t! ^ 


o g 

1 ^ ^ 

cd ^ £-j 
”■2 3 

2 .2 HH 

^ S ^ 

& bOxJ 

§ 12 o 


^ £ I 

^ ° I 

.3’2 O 
^ c3 ‘a 

T3 cd 

a 00 'P 


^ <u *-* 

S o ^ 

' O d ^ 

6b O 

(D 


^ s: ^ 2^n 
S o 


^ U ^ c/3 ^ 

|l ^ I g^ 

£ 1.2 
^ o> a ^ . 
a d o 
o ^ o o ^ 
•2 S3 ^ ^ 


.SP-^ § S P <" 

w O c2 cd J3 o 


Knowledge of soil properties and below ground C stores is limited. While optical sensors do not penetrate 
into the soil, some surface properties can be estimated. Ahn et al. (1999) used AVIRIS data to map soil surface 
properties, demonstrating the potential to map heterogeneity and identify localized characteristics, e.g., soil organic 
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Asner, G.P., C.A. Wessman, and D.S. Schimel, Heterogeneity of savanna canopy structure and function from 
aging spectrometry and inverse modeling. Ecological Applications 8:926-941, 1998. 
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application in crop growth models for predicting leaf senescence. Annals of Botany 85:579-585, 
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MAPPING ARID VEGETATION SPECIES DISTRIBUTIONS IN THE WHITE MOUNTAINS, EASTERN 
CALIFORNIA, USING AVIRIS, TOPOGRAPHY, AND GEOLOGY 
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canonical space were calculated. Results were expressed in standard deviation units - closer to the centroid means 
higher probability of occurrence, and beyond three standard deviations was taken as absence. Axes maps show the 
spatial realization of the linear combinations of environmental variables. Each cell has six axis values that are then 
used to determine whether that cell falls into environmental envelopes of selected species based on the species’ 
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Axes and species parameters were derived from 202 calibration plots. 146 plots were used to evaluate 
presence/absence of 17 species. Projections were smoothed by 3x3 cell moving average (to account for registration 
errors between the different environmental and spectral grids). Kappa statistics were calculated for each species. 
Kappa is the relative increase in correct assignments over that expected by chance. A threshold standard deviation 
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Menodora spinescens Pinus monophylia 




t: tt o 

c /3 Oj ^ 
T3 O 

n 8 g 

.2 ft ® 

^ W c /3 
W CO 

^ !Hj 

C/3 

« .3^ 

If ^ 

' ^ cn 

^CT) ^ 


O .1=5 . 

^ P 'Z:' 

gl a 

.2 <u s 

§> 2C 

•S Jh W 
<D .1=^ 

a s 

K (D 0) 

i| ^ 

Cl 

^ rr» ^ 



(D I" I 

+-» ca o 

5> 13 

^ t/3 

^ S 3 

« ” o 

■S § ^ 


}rJ <i> 

I s 

•Jag 

h U 22 


o « £P 

III 

ll^ 

^ n-i 

•o .2 ^ 

S Qi ^ 

'S O S gB 

al ^ ’R 
.. § S S. 


^ o 
R Q S 
a« S ^ > 


i |i 

/--N H <L> 

a? -5 bo 

•S s s 

Ph M 2 

p ^ 

1^ I 

c c /3 o 
.P (D ^ 

fe ^ trt 

'w' Q 03 

s§ S 

^ W) Ph 
^.S ^ 

^ o "w 

s '? g> 

0 ji .S 

^ !=^ tP? 

S 00 ± 

,R czi (li 

‘ft §^ ® i 

O g 7d ^ 

0 0^1 

1 's' > 

■i ci^ o 

S 6 ^ 


a« *S ^ 

® S'O 2 

^ S go 

2 ■§ -ft 

® i I 8 

.„sfli i 


S o ^ 

^ S O O 
O Co cn (D 

i 

a cT s ^ 

ft § g 

S « 1 J 

.2f B ^ p 

^ P <D 

o o :g S 
H S S 


ft a 2 
^a g 2 - 

I O c 3 
<1^ rn 

•S s ^ 


B 5P^ 

I O 
TJ ^ VO 
1) V-i CO 
^ O w 

S ft" 
i =* g 

a <D .S 
O U gr, 

O iJ W3 

«3 X i 

.^ <D Vh 

^ in 

> *z: ^ 

O o3 o 

^ O 1 
O ft 2 

S- ^ ® 

. -S ^ 

*3 O bD 
R tP? P 
O '*7' -S 
s? ft ^ 

3 '? 

^ o ii 
R i/j 

^ C« S 

cP g 03 Qj 

’ft R Q J 

g I ® "O 
•2 2 .§ 

.§ §§1 
B ^ CO C^ 

.52 R ^ 

^ ^o P 
o s*-H g -b 

•»-! > c/a 

'P . (D o-j 

a S V s 

ft|^ =« 

a I ■§ I 

p -a 2 

*9 ^ w 

9 '*■’ CO 4 ^ 

o p tj ^ 

JX rP n 

W 03 03 ci= I 


§ « O o 
§ 8 
^^g I 

R fo S 

so .« to 

I °- § ■§ 

%lll 

p 2 S ^ 

■3 a V2 

^ 3 ^ " 

.52 O c3 P 

'O g CO O 

(D o ^ 

^ ■s g. § 

^ a « ‘S 

CD O .— • CO 


s •§ a 

WD O o 

CO VO 

9 § 









3. EVALUATION AND RESULTS 
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factors was effective in determining the fundamental temperature/moisture gradient. 
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ANALYSIS OF SUB-PIXEL MIXING IN HIGH-ALTITUDE AVIRIS DATA OVER 
VIRGINIA CITY, NEVADA, USING SYSTEMATIC FIELD-BASED OBSERVATIONS 
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3. Field Methods 



Figure 2. Comparison of spatial resolutions for AVIRiS (18-m) and DOQ f1-m) The SW 
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complementary data sets for the unambiguous identification of sulfates, clays and silicates. 
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rsummS“n 9^^^^ Imaging Spectrometer (AVIRIS) data were collected over the site 

le summer or i i his type of remotely sensed data has a snatijil ^ 
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Figure 2~Generalized geologic map of the Hillsboro mining district (McLemore et ah, 2000). 
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A model for hydrothermal alteration around porphyry copper deposits was developed by Lowell and 
Guilbert (1970). In this model, alteration progresses outward from a central core of potassic alteration to a phyllic 
alteration zone, then to an argillic zone, and finally to a propylitic alteration zone along the periphery. Potassic 
alteration is characterized by biotite and orthoclase. Phyllic alteration consists of quartz, sericite and pyrite. Argillic 




alteration mineral include kaolmite and montmorillonite. Typical propylitic alteration minerals include epidote 
calcite, and chlorite. A schematic cross-section through an idealized copper deposit is illustrated in Figure 3. 
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Figure 4— Preliminary map of selected minerals in the Hillsboro minin g district. 
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SPECTRAL BAND CHARACTERIZATION FOR HYPERSPECTRAL MONITORING OF WATER 

QUALITY 

Stephanie C. Vermillion,' Rolando Raqueho,^ Rulon Simmons^ 
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absorption curve as a proportional indicator, and all other concentrations were scaled to that curve. 
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3. EXPERIMENTAL RESULTS 
3.1 Band Averaging 




carbon (CDOM), and 2.32 mg/m^ (CHL). TSS was the least affected by changes in spectral resolution, ’having the 
smallest error between resolutions ( x = 0.26) as well as the smallest range of error (1.07). The error ranges for 
CDOM and CHL were 2.19 and 2.52 respectively. 
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Table 4. Band loading for the first principle component of the 30 
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Comparison of AMOEBA Predicted Concentrations vs. 
Ground Truth Measurements 
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The last step prior to analysis was to add random noise to the dataset to simulate various SNR levels. This 
step was performed using equation 1 from Okin and Roberts (2000): 



2.3 Simulated Product Assessment 
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Figure 2. NDVI images with varying spatial resolution, spectral resolution, and noise levels. Visually, spatial 
resolution is more important than SNR when SNR for NDVI bands is greater than or equal to 100:1. 
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CORRELATING AVIRIS IMAGERY TO FIELD SAMPLING AND SPECTROMETER 
MEASUREMENTS FOR INORGANIC SOIL CARBON 
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/'.I -1 ® ^ nave ucmonsiratea me abilitv to accnratplv pctim 

many of the soil p^ameters, such as salinity, texture, organic matter and associated compraents (Ben-DOT and 
Bamn 1994 Condtt, 1972; Csillag et ah, 1993; Stoner and Baumgardner, 1981; StonerTt ari98m ^ 
important studies have fr^sferred these techniques of soil characteristic identification to airborne andlateL 
images (Mettemicht and Zinck, 1996; Palacios-Orueta and Ustin, 1996). 



varying salt concentrations spatially. 
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Presented here is a preliminary review of the soil spectra, chemistry and mineralogy in this study. The 
Division of Natural Resources and Agriculture laboratory, UC Davis, laboratory conducted the chemical analyses on 
the 63 soil samples. The laboratory spectra was collected using a LabSphere diffusion sphere attachment to a Cary 5 
laboratory spectrophotometer (Varian Inc.), at the Center for Spatial Technology and Remote Sensing, UC Davis 
(CSTARS). The soil samples were prepared by lightly grinding in a mortal and pestle, then removing remaining 
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parameters and viewing geometry values were measured in the field at the time of data collection for 
Pinus strobus: tree density, 1100 trees/ha with Poisson distribution; irregular ellipsoidal crown shape, 
height of trunk, 8.5 m; height of tree, 15.0 m; trunk radius, 8.3 cm; crown radius, 2.0 m; spherical leaf 
angle distribution; shoot area, 0.0008 vci\ canopy effective LAI, 2.5; and leaf area density, 0.4171 /m. 



demonstrates that radiative transfer model inversion techniques can be applied successfully to 
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Canopies. Remote Sensing of Environment. BOREAS Special Issue III (submitted, 2001b). 



